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FOREWORD 


This report reveals a feature of the noise generation process in turbofans described herein as 
"mode trapping". As with most physical phenomena, its understanding developed on more than 
one front and several people contributed to the evolution of the idea. Test data characteristic of 
mode trapping and elements of its physical explanation were first identified in a paper by Topol, 
Holhubner, and Mathews entitled "A Reflection Mechanism for Aft Fan Tone Noise" (AIAA 
Paper No. 87-2699, October 1987). At that time, the author of this report (Hanson) was 
employed at the Hamilton Standard division of United Technologies and working on noise 
associated with unsteady coupling between rotors of a counter rotation prop-fan. While 
consulting to Hamilton Standard, also in 1987, Professor Marten Landahl of MIT recognized the 
possibility of trapping a wave between 2 rotors and proposed a method for computing the effect 
via an impedance model. In conjunction with Landahl’s work. Professor Kenneth Hall, then at 
UTRC and now at Duke University, originated the idea of using the S. N. Smith’s unsteady 
cascade code to represent the coupling and developed a model for the interaction of 2 unloaded 
rotors at a single harmonic. 

In 1990, at the request of the Pratt and Whitney division of United Technologies, Hanson 
continued work on the reflection mechanism to put it on a firmer analytical basis for ducted fans. 
His work was funded by NAS A- Lewis in 1991 as part of Contract NAS3-25952, during the 
course of which the author transferred to Pratt and Whitney. NASA project management was 
provided originally by Dr. Christopher Miller and finally by Mr. Dennis Huff. The author owes 
thanks to all of the individuals named above for their support and technical assistance. Thanks 
are also due to Professor Edward Kerschen of the University of Arizona who provided valuable 
encouragement on use of the actuator disk model employed herein and on the matrix method 
finally used to solve the coupled equations. 
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SUMMARY 


Typical analytical models for interaction between rotor and stator in a turbofan analyze the 
effect of wakes from the rotor impinging on the stator, producing unsteady loading, and thereby 
generating noise. Reflection/transmission characteristics of the rotor are sometimes added in a 
separate calculation. In those models, there is a one-to-one relationship between wake harmonics 
and noise harmonics; that is, the BPF (blade passing frequency) wake harmonic causes only the 
BPF noise harmonic, etc. This report presents a more complete model in which flow tangency 
boundary conditions are satisfied on 2 cascades in relative motion for several harmonics 
simultaneously. By an extension of S. N. Smith’s code for 2D flat plate cascades, the noise 
generation/frequency scattering/blade row reflection problem is solved in a single matrix 
inversion. It is found that the BPF harmonic excitation of the stator scatters considerable energy 
into the higher BPF harmonics due to relative motion between the blade rows. Furthermore, 
when swirl between the rotor and stator is modeled, a "mode trapping" effect occurs which 
explains observations on fans operating at rotational speeds below BPF cuton: the BPF mode 
amplifies between blade rows by multiple reflections but cannot escape to the inlet and exit ducts. 
Flowever, energy scattered into higher harmonics does propagate and dominates the spectrum at 
2 and 3 times BPF. 

This report presents the complete derivation of the theory, comparison with a previous 
(more limited) coupled rotor/stator interaction theory due to Kaji and Okazaki, exploration of the 
mode trapping phenomenon, and parametric studies showing the effects of vane/blade ratio and 
rotor/stator interaction. For generality, the analysis applies to stages where the rotor is either 
upstream or downstream of the stator and to counter rotation stages. The theory has been coded 
in a FORTRAN program called CUP2D, documented in Volume II of this report. 

It is concluded that the new features of this analysis - unsteady coupling, frequency 
scattering, and flow turning between rotor and stator - have a profound effect on noise generation 
caused by rotor/stator interaction. Treating rotors and stators as isolated cascades is not adequate 
for noise analysis and prediction. 



SECTION 1 
INTRODUCTION 


Rotor/stator interaction is undoubtedly the oldest and most frequently treated subject in 
turbofan noise research. Countless reports and journal articles have appeared over the past three 
decades analyzing the complex aerodynamic interference and acoustic reflection effects that 
contribute to noise at blade passing frequency and its harmonics. Geometry of interest shown 
in figure 1 is representative of modem engines where the rotor is upstream of the stator. The 
first major insight into the acoustics of this kind of fan was provided by Tyler and Sofrin (ref.l) 
and then Sofrin and McCann (ref. 2) who revealed the link from number of blades, number of 
vanes, and tip speed to the duct modes and their propagation/decay behavior. Methods to 
compute the strength of these modes via unsteady aerodynamic calculations at first used crude 
representations of unsteady loading response at the stators. Then Ventres, Theobald, and Mark 
(ref. 3) at Bolt Beranek and Newman developed a computer code (herein called the BBN code) 
that linked a 2D flat plate cascade theory for the unsteady aerodynamics to the 3D annular duct 
acoustic modes via a strip analysis. Their theory included axial flow and non-compact treatment 
of a stator that could be leaned, tapered, or swept. The unsteady aerodynamic model in the BBN 
code can be explained with reference to figure 2. Rotor blades trail viscous wakes producing an 
upwash distribution on the stator. The stator responds with the lift response calculated to produce 
an equal and opposite upwash so as to satisfy the flow tangency boundary conditions on the 
vanes. The computed loading then acts as a dipole distribution in the ensuing coupling to the 
duct modes. In figure 2, the rotor blades are shown with dashes to emphasize that their role is 
only to generate the wakes; any reflection effect that they might have on the acoustics or 
unsteady aerodynamics is ignored in the BBN code. 

Nevertheless, the BBN model includes much of the important physics of the problem and 
has been exercised extensively by Topol (ref. 4) in comparison with test data. However, in 
another study, Topol, Holhubner, and Mathews (ref. 5) found that the method failed to predict 
observed phenomena, even qualitatively at some RPM’s, for a fan designed for BPF cutoff. 
These authors noted that BPF modes could cut on in the swirl !“ >w region between the rotor and 
stator at speeds considerably lower than in the inlet and exit flow. In this intermediate speed 
range (between BPF cut-on in the swirl region and cuton in the aft duct flow), 2xBPF and 3xBPF 
harmonics (which are cut on) in the aft region were found to increase to levels far above those 
predicted by the BBN code. Aft-radiated sound power behavior for the first 3 BPF harmonics 
is shown in figure 3 for a typical, current day, high bypass ratio fan engine with 38 blades and 
72 vanes. The revealing behavior to note in figure 3 is the rapid onset of 2xBPF and 3xBPF at 
the RPM where BPF cuts on in the swirl region. Reference 5 concluded that the increase in 
higher harmonics in the aft region was due to the following sequence of events. The rotor wake 
excites the stator predominately at BPF; this sends a BPF wave to the rotor that scatters it into 
higher harmonics; these harmonics are cut on and propagate through the stator to the far field. 
This hypothesis was supported by analytical modeling that included various simple models for 
the source strength and for reflection and transmission coefficients at the rotor and stator. 

The BBN code cannot be used to check this reflection/scattering hypothesis because it 
includes neither the frequency scattering mechanism nor reflections from the rotor. In the present 
paper, the model depicted in figure 2 is refined to include the rotor in the aerodynamic 
calculations via a fully coupled unsteady analysis satisfying the How tangency boundary condition 
on both blade rows simultaneously for several harmonics. Thus, noise generation, frequency 
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scattering, and reflections from both rotor and stator are solved simultaneously via a single matrix 
inversion. Furthermore, swirl between the rotor and stator is modeled via a simplistic system 
that turns the flow at the rotor leading edge and straightens it at the stator trailing edge. (This 
is a first approximation to including blade mean loading in the analysis by concentrating it at 
actuator disks.) Modal jump conditions across these disks are derived to satisfy continuity of 
mass and momentum. It should be pointed out that Kaji and Okazaki developed a coupled 
rotor/stator method (ref. 6) in 1970; however, this model was confined to one frequency and did 
not include the swirl effect. Correlation with their work is shown later in this paper. 

This work is an extension of the Smith code (ref. 7), which is an industry-standard 2D, 
compressible, linearized, unsteady flat plate cascade method. The coupled cascade theory has 
been coded in FORTRAN 77 as CUP2D and developed on UNIX™ workstations. The theory 
applies to counter-rotation stages and to rotor/stator stages with the rotor either upstream or 
downstream of the stator. Applications herein are primarily to rotor/EGV combinations, but 
comparisons are given with the Kaji and Okazaki results for IGV/rotor interaction. 
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SECTION 2 

REVIEW OF SMITH’S CASCADE METHOD 


Since this paper is an extension of Smith’s theory (ref. 7), we start with an overview of his 
method. Further description can be found in a survey paper by Whitehead (ref. 8) and, of course, 
in the original Smith paper. Figure 4 shows the basic geometry to be a single cascade of flat 
plates aligned with the mean flow direction. The objective is to find the unsteady loading 
distribution on the cascade (and the resulting output waves) for a specified velocity disturbance 
input to the system. Smith’s approach is to derive an integral equation for the upwash on the 
cascade caused by the unknown load distribution, discretize the equation, and solve for the loads 
by matrix inversion. 

Smith’s analysis starts with the linearized differential equations for continuity and 
momentum. He finds solutions for velocity and pressure disturbances in the form 

u 

- e i(ax+Py +co/) (1) 

P. 

where the space, time behavior is given by the exponential and the overbarred quantities are 
complex amplitude coefficients. The kernel of the integral equation is expressed as a Fourier 
series in the y coordinate with terms counted by an index k ( r in Smith’s notation), which 
corresponds to the circumferential spinning mode index in the Tyler-Sofrin theory. Transverse 
wavenumber (5 can take on all values given by 

(5 = ^ IlHlIL , -°°<k<+°° ( 2 ) 

where c is the interblade phase angle which will be discussed later, a and frequency to are 
inputs to the problem. 

Substitution of equation 1 into the original continuity and momentum equations leads to a 
characteristic equation with two families of solutions, each having a different form for the axial 
wavenumbers. In one family are the vorticity waves which are carried at the mean flow speed 
in the downstream direction only. Components of these waves are denoted by the subscript 3 
as in u 3 , v 3 , and p 3 . Axial wavenumbers for the vorticity waves are 

= ( 3 ) 

3 U 

The second family are the pressure waves, which propagate at the speed of sound. Upstream- 
and downstream-going pressure waves are denoted by the subscripts 1 and 2, respectively and 
their axial wavenumbers are given by 

(4) 

where a is the speed of sound and the square root contains the cutoff information critical in fan 
acoustics. When the argument of the square root is positive, the axial wavenumber is real and 


_ t/(0) + Vp)±a\/(0) + Vp) 2 -{a 2 -U 2 )$ 7 


a 2 -U 2 
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the wave propagates at constant amplitude in the duct; in this case, the ± choice is made to give 
a shorter wavelength upstream and a longer wavelength downstream. When the argument of the 
square root is negative, the wave decays exponentially; the ± choice is made to guarantee 
vanishing disturbance for x—*». The cuton criterion can be simply stated: a wave with 
transverse wavenumber p will propagate if the frequency is "large enough", that is if (co+Vp) 2 
>(a 2 -U 2 )$ 2 . The higher the mode order, the higher the frequency required for propagation. 


Blade loading is in terms of the vorticity distribution T(z 0 ) over the source coordinate 
z 0 (along the chord). With this notation. Smith wrote the following for the transverse velocity 
in the pressure wave anywhere upstream of the cascade 



E 


k 


/ i(a,Jt+py + 0 )/) -f(a,cos0 +Bsin0)z„ , 

Vj e 1 e ° dz 0 


(5) 


where c is the airfoil chord, v/ is a function of co, P, 0, and W given by Smith and defined 
in Appendix C. Expressions for the downstream velocity components are obtained by replacing 
the 1 subscripts with 2’s and 3’s. Axial velocity and pressure are linked to the v'’s as follows 


/ 

“l,2 



( 6 ) 


and 


P\,2 = -(C0*t/«| 2 + yp)p„ 


(7) 


for pressure waves and 



(8) 


for the vorticity waves (which do not produce a pressure disturbance). 


Finally, to discretize equation 5, the load distribution is divided into N p panels and the 
integral is converted to a sum according to 




N„ 


The result is 


N p OO 

W v-' f /(a.jr+P y+M) -i(a,cos0+psin0)z • 

Vi= — E E V 1 e e JxL( j > 


(9) 


(10) 


sic J _| £=-00 

where z 0 j is the normalized chordwise coordinate. 


The next section shows in general terms how this type of expression can be adapted to the 
blade row interaction problem. 
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SECTION 3 

DEVELOPMENT OF THEORY FOR LOADING AND ACOUSTICS 
OF COUPLED CASCADES 


The previous section described the mechanics of finding waves generated by the unsteady 
loading distribution on an isolated, infinite (in the y direction) 2-dimensional cascade. Geometry 
to be analyzed for the coupled cascade problem is sketched in figure 5. The analysis is kept 
general so that either blade row can rotate in either direction; thus, counter-rotation cases and 
IGV/rotor combinations are included. In spite of this generality, the upstream blade row will be 
called the rotor and the downstream blade row will be called the "stator" to simplify the 
terminology in the lengthy derivation that follows. The typical case of interest herein is for an 
upstream rotor and a downstream stator. In this case, the downstream blade row would be fixed 
and 0j would have a negative value. If the downstream blade row is moving, its tangential 
speed is input as a negative number. 


3.1 DUCT COORDINATES AND MULTIPLE HARMONICS 

The 2-dimensional geometry of figure 5 will be adapted to annular geometry (annulus radius R ) 
for the case of a "stator" with B 2 rotating at angular speed Q 2 ar| d a rotor wit h B i blades 
rotating at angular speed £2j. The stator is fixed in the x,y 2 coordinates with inflow angle 
0 2 defined by its rotational speed £l 2 R and the swirl velocity V y as shown by the velocity 
triangle. The rotor is fixed in its x,y l system which is moving up at speed Q.]R. The relative 
air speed in the tangential direction is Q.jR-V s which defines the inflow angle 8y. 

In this section we adapt Smith’s formulas to the notation of duct acoustics and account for a 
general periodic waveform rather than a single harmonic. This is done by examination of the 
term in equation 10 that we call the kinematic phase 

\|/ = (3y+cor 

which is discussed in some detail in appendix A. Frequencies of stator excitation, of course, are 
the blade passing harmonics, co = nBjQ,, where we have defined the relative angular velocity 
£2 =Qj-Q 2 . We normalize the transverse coordinate _y 2 by the radius of the annulus R, 
effectively converting to cylindrical coordinates with y 2 = R( j> 2 for the stator. Appendix A 
shows that the normalized transverse wavenumber can then be written as 

PR=P„ t =-(*B, -*B 2 ) 02) 

so that for the stator 

= (13) 

This is the exponential already familiar from the Tyler-Sofrin spinning mode theory (ref. 1). 
Equations 12 and 13 represent a mode with nB I -kB 2 lobes spinning at angular speed 
nBjQ/fnBj-kB?), as viewed in the stator frame of reference. 

Given this notation, equation 10 can be generalized immediately for stator waves with 
multiple harmonics: 
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V 1 = M r 2 £ £ £ v{ * ^ , - 1 (a ^ COsQ l x L(nj - } (14) 

a SC 2 „=i /= i it = _oo 

Overbars denote normalization by /?, the dimensionless axial wavenumber is a y /? = a lnk , M r2 
is the relative Mach number W^a, and SC 2 is the spacing to chord ratio for the stator. 
Working forms for the \|/’s are given in appendix A, for the a’s in Appendix B, and for the 
v'’s in appendix C. 

As a preview of the frequency scattering phenomenon, note that y nk for the stator waves 
can be viewed in rotor frame of reference via the transformation <J> 2 =<f> y +£2f. The result is 

Vnik“M i*kB 2 at (15) 

This shows that the waves caused by any one stator loading harmonic n appear at all multiples 
of vane passing frequency in the rotor coordinate system. 


3.2 ACCOUNTING FOR FLOW TURNING AT THE ROTOR AND STATOR 

The model discussed so far includes swirl in the flow entering the stator but no turning of 
the flow by the rotor and stator loading. Since swirl was fundamental to the mechanism 
proposed by Topol et al (ref. 5), a method was developed to represent turning of the flow that 
could be combined easily with elements of the Smith code. The scheme is to use actuator disks 
to turn the mean flow at the rotor leading edge and at the stator trailing edge as sketched at the 
top of figure 6. Since the mean flow properties (tangential and axial velocity, density, speed of 
sound) jump across the inlet and exit interfaces, the unsteady quantities jump also. This is 
handled by deriving reflection and transmission coefficients for each type of wave to maintain 
continuity of mass and momentum (linearized) on a mode-by-mode basis. Derivation of these 
coefficients is given in appendix D. 

Since the actuator disk representation of the mean loading effect seems to be new in this 
context, a few more words are in order. We consider the blade loading to be the sum of a steady 
part and an unsteady part. The unsteady loading is distributed over the blade chords via panels 
as in the original Smith method and is considered small enough to be treated as a perturbation 
of the mean flow. The steady loading is concentrated at the rotor leading edge and at the stator 
trailing edge and may be large enough to have a significant influence on the mean flow. In fact, 
this is its main role in the analysis; it provides the background flow that is perturbed by the 
unsteady loading and controls the propagation properties for the various types of waves. Thus, 
we can think of building up the loading as follows. First, turn on the mean loading. The mean 
flow enters the stage axially, turns at the rotor leading edge, and straightens out again at the 
stator trailing edge. The mean flow is thus divided into three regions by the two actuator disks. 
Mean flow properties in the three regions are considered known from a separate analysis. Next, 
the unsteady loading is activated; the waves behave according to Smith’s theory in each region 
with jumps across the actuator disks to preserve continuity and momentum. 

Smith’s analysis was based on the assumption of sources generating waves in an infinite 
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uniform medium so that there were only outgoing waves. These waves were solutions to a set 
of inhomogeneous equations with blade vorticity representing the loading source. Now we must 
include waves that are solutions to the homogeneous equations to account for the reflections from 
the inlet and exit interfaces. The situation is sketched in figure 6 where we show the loading 
source LS(n,j) sending out the direct source waves KS1, KS2, and KS3 from Smith’s analysis. 
The reflected waves VS1, VS2, and VS3 will be derived in the next section in terms of the inlet 
and exit reflection coefficients. With the KS's and VS's as known functions proportional to 
the loading source, the entire wave field is still driven by the unknown loading elements LS(n,j) 
as with the original Smith method. 

For the reflection and transmission coefficients, three interactions had to be considered, as 
shown in figure 7: upstream-going pressure waves Vj could transmit a pressure wave v 4 and 
reflect a pressure wave v 2 and a vorticity wave v 3 . Downstream-going pressure v 2 and 
vorticity v 3 waves could each transmit pressure and vorticity waves (v 8 and v 9 ) and reflect 
pressure waves (v;). It was most convenient to derive the coefficients in terms of the modal 
transverse velocity components as follows. 

Upstream-going pressure waves at inlet 

v_2 la bhk = R l lW)V\ la b$nk 

v 2 ,a b a ink =R l^ n ^ v _l /a b^nk (16) 

v 4 /a aKk = T l4 n ^ v l ,a bhk 


Downstream-going pressure waves at exit 

7, * *"**•_ la b fa = R 2l (n,k) v 2 e ia ™?'la b fa 
~8 * ^yaAk = T 2 ,(n,k) v 2 e ^ !a b fa 
~9 e lCL9nkXe/a c a 9nk = T 29^~2 eia2 " kXe,a b^nk 

Vorticity waves at exit 

v, « ia '" kX y$nt - *3, W)Vj« !a b<hnk 

v 9 e ia9 " kXe l«c a 9nk = T 39 (n,k)v 3 e lCLsnkXe !a b a 3nk 


(17) 


(18) 


Smith’s subscripting convention has been continued so that 4 denotes the pressure wave upstream 
of the rotor and 8 and 9 denote the pressure and vorticity waves downstream of the stator. 
Derivation of the reflection and transmission coefficients is given in appendix D. It will be shown 
later that reflections are small at interfaces between regions where a mode is cut on. However, 
for modes passing from a cut on region to a cut off region, the reflection coefficient is of order 
unity. Thus, modes can become trapped between the rotor and stator and can be amplified, as 
will be seen. 
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3.3 EQUATIONS FOR COUPLING 


The wave field between the rotor leading edge and the stator trailing edge can now be 
constructed in terms of direct waves from the load elements (via Smith’s methodology) plus the 
waves reflected from the inlet and exit interfaces as sketched in figure 6. Equations coupling the 
rotor and stator will be derived in 4 arrays of influence coefficients as follows. 

KSS - Stator on Stator 
KSR - Stator on Rotor 
KRR - Rotor on Rotor 
KRS - Rotor on Stator 

Later, these will be assembled into a large matrix to represent the entire coupled system and the 
system of equations will be solved by inverting the large matrix. The difficult parts of KRR and 
KSS have been taken care of by existing subroutines from Smith’s original code. In particular, 
all issues of singularities that occur as source and field points approach each other are handled 
by these routines. KRS and KSR and the parts of KRR and KSS associated with the reflected 
waves are new and their derivation is described below. 


Waves Caused by Stator Unsteady Loading 

Consider first the direct waves caused by the stator loading. In terms of a kernel KS1 
defined below, we can write for the upstream wave of equation 14 

v ldlr *<» t E E KSl(n,kj)xLS(nj) (19) 

n- 1 lc=-» j= 1 

Similarly, the downstream waves are 

v 2dir ~ a b E E V nlt e* a ^ x ~^'^'E,KS2(nJ<j)xLS(nj) (20) 

n - 1 k - -°o j-\ 

OO OO — — N p 

v 3dlr ="(. E E KS3(n,kj)xlXnj) (21) 

« = 1 &=-» j - 1 


where the kernels are defined by 
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M r2 v{ 

KS1 = _ - J- e 

SC 2 fink 

M r2 v 2 

KS2 = — — — — e 

fink 


-i(a Ink cosQ 2 +V nk sinQ 2 )c 2 z OJ 
-i(a 2nk co$Q 2 + P„*sin Q 2 )c 2 z oj 


KS3 ~ 


M 


r2 v 3 - '( a*,*cos 0 2 + P„* sin 0 2 ) c 2 z oj 


SC 2 


( 22 ) 


(Use of v[l§ nk , v[lfi nk , and v 3 / /a 3 ^ avoids problems later with divisions by 0.) The 

reflected waves are written in parallel form in terms of the coefficients VSl, VS2 , and VS 3 which 
will be found in terms of the reflection coefficients derived in appendix D: 


oo 

oo 


N P 


v lref = a b E 

E 


*(«/„** + V„*)^ VSl(n,kJ)xLS(nj) 

(23) 

/! = 1 

k- 




oo 

oo 


N p 


v 2ref = a b E 

E 

P/ifc e 

' (a 2^ + ^E VS2(n,kJ)xLS(nJ) 

(24) 

n - 1 

/C=-oo 


M 


oo 

oo 


n p 


v 3ref ~ a b E 

rt = 1 

E 

/c= -oo 

a ?,nk ' 

e ***** +%,k> E VS3(n,kJ)xLS(nJ) 

(25) 


Now, by applying the definitions of the reflection coefficients in equations 16-18, 3 
equations can be written for the 3 VS coefficients: 

VSl e iainkXe = ( KS2 e + VS2 ) /“ w " s /? 21 W) 

+ (KS3 e " iaiHkXs + VS3 ) e R 3l (n,k) (26) 

VS2 = (KS1 e ~ ia,nkX * + VSl ) R l2 (n,k) 

VS3 = (KS1 e~ i<X,nkXs + VSl)R u (n,k) 

This set of equations can be solved for VSl, VS2, and VS 3 in terms of KS1, KS2, and KS3 and 
the reflection coefficients, effectively determining the entire wavefield generated by any stator 
loading element LS(n,j). The results upstream and downstream of the stator are 


up 


= a bY, E E e Nnk [$ nk {KSle ia,dX - Xs) + VSle ia,nkX + VS2e l(X2nkX ) 


n=\ k= -<» 7=1 


(27) 


+a 7ink VS3e iaMkX ]xLS(nJ) 
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and 


v dn= a bT, E E ‘ ‘ V “ 1 P„*(SS2<! * VS I e ia ”*' * VS2e ia ’' tX ) 

rt=l *« — y=l ^ 

+ a 3nk (KS3 e ~ Xs) + V53 « ] x L5(n j) 

These velocity fields have only to be evaluated at control points on the stator and rotor to find 
the effect of the stator on itself, KSS, and the effect of the stator on the rotor, KSR. This is 
treated in the next two sections. 


Effect of Stator on Stator 

We place one control point on each stator loading panel to give the values for x and <J> 2 
for use in equations 27 and 28. Control points, counted by i from 1 to N p , are at the same 
chordwise locations z t - used by Smith and lead to the following: 

X-*X s +C 2 Zi ;COSe 2 (29) 

<|> 2 — » c 2 z,sin0 2 


In examining equations 27 and 28, note that the terms with the KS’s are the direct waves that 
Smith already accounts for in his method. The remaining terms, with the V5’s are common to 
both equations and can be evaluated at the control points. To find the upwash (normal to the 
stator chord), we add the transverse component times cos0 2 and minus the axial component times 
sin0 2 , the latter requiring Equations 6 and 8. For control point i the result can be written 


W 2 


= E '£KSS , (n,iJ)xLS(nJ)e inB,il ‘ 

n-\ y=l 


(30) 


where the primes on w ssi and KSS are a reminder that this part of the stator upwash is due to 
the reflected waves and is to be added to matrices from the original Smith method. The array 
of influence functions is 


oo 

KSS'(n.ij) = T J_x 

Kp„ t cos8 2 -o^sinBj) VS l e ' (a “ cos9 ^ |, ”* slne ^ z ' 

+(P„ t cose 2 -arsine,) VS 2 e‘ a ^. e ^«*<>i 
+(a 2nfc cos6 2 + |5„ 4 sin0 2 ) VS3e ta ^‘e , 


Equation 30 gives the waveform of the upwash at the stator control points z- t as an infinite 
Fourier series. For a matrix solution, we can satisfy the flow tangency requirement 
simultaneously for only a finite number of harmonics, N h . (By adopting the Smith code, we 
have already accepted tangency at only a finite number of control points, N p , along the chord.) 
Since each of the loading harmonics on each blade row induces upwash at the other blade row 
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in all harmonics, it is not clear a priori that truncating the harmonic series will lead to acceptable 
results. It does turn out, however, that the computed results for the lower harmonics are not 
strongly influenced by higher harmonics. A demonstration of this is given later in Section 4. 
Based on the above argument, the upwash waveform can be written 


w 


/ 

ssi 


W 2 



WSS'(n,i) 


inBiLIt 

e 


(32) 


where 

n p 

WSS ' (i n,i ) = £ KSS ' ( n,iJ)*LS(nj ) < 33 > 

7=1 

When KSS / is added to the matrix from the direct waves computed with the Smith subroutines, 
that matrix is called KSS and the net effect of the stator loading elements at the stator control 
points is 

N P 

WSS(n,i ) = £ KSS(n,iJ) x LS(nJ) ( 34 ) 

7=1 

This is the first of the coupling equations. It can be seen that, for the effect of the stator on 
itself, each loading harmonic n couples only to the same harmonic in the upwash. 


Effect of Stator on Rotor 

For this effect, we proceed in a similar fashion. But in this case, we will find that each 
loading harmonic on the stator couples to all of the upwash harmonics on the rotor. In the rotor 
coordinate system the control points are 

x = c 1 z / cos0 ] 

<!>! =Cj z ( sin0 1 


which can be applied to the form for V|/ nfc given by equation 15 or A- 17. These apply to waves 
caused by stator loads as viewed in the rotor system and explicitly shows the scattering from 
loading harmonic n into the other harmonics counted by k. The result, written for a finite 
number of harmonics, is 


w Nh Nh N p , D ^ 

E E KSR(n,Uj)xmnj)e' kB!a ‘ 

W l n = 1 k=-N h ;= 1 


(36) 


where 


KSR(n,k,ij ) = 

—L[ (P^cose, -a lrt sin9,)(ra/ +KSle 
M r i 

.(P„ t cos6, -arsine,) VS2 e 

♦(choose, *p„,sine,) vw/ (a -« cos9 '‘ |) "* sin9 ' )r ' z ' 1 
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Note in equation 36 that negative frequencies have appeared as a result of viewing stator waves 
in the rotor coordinate system. To obtain a form parallel to equation 30 with only positive 
frequencies, we rewrite equation 36 using the following strategy. The ^-summation can be 
written as a single sided series: 


w 


srt 


W, 



£ 52 [KSR(nAiJ)xLS(nj) e + KSR(n,-k,iJ)xLS(nJ) e ~ lkB 2 at ] (38) 

n=l ;'=1 


But only the real part of this expression has meaning so we are free to take the complex 
conjugate of any term. Applying this to the second term above and denoting conjugate by a *, 
we find 

N h N h N P 

— " E £ £ { KSR(n,k,ij) x LS(nj) + [KSR{n,-k,ij)xLS(nJ) ]* (39) 

w i *= l «= l j=i 

which is in terms of positive frequencies, as desired. Now, in parallel with Equation 30 for the 
stator upwash, we can write the upwash waveform as 

— = 52 WSR(k,i) e ( 4 °) 

w i k* l 

where the upwash harmonic on the t th rotor panel is given by 


N h "p 

WSR(k,i) = 52 52 { KSR(n,k,iJ) x LS(nj) + [ KSR(n, - k,ij ) x LS(nJ) ] * } < 41) 

n=l j= 1 

This is the second set of coupling equations. Again, note that each stator loading harmonic n 
contributes to all upwash harmonics k on the rotor. 


Waves Caused by Rotor Unsteady Loading 

The derivation for the effect of the rotor loading on the rotor and stator proceeds along 
similar lines; the procedure is only outlined here but the working equations are presented in 
detail. The main differences between the formulas for rotor and stator loading waves is 
permutation of the 1 and 2 indices, permutation of the n's and &’s, exchange of /?’ s and S ' s in 
some of the symbols, and the appearance of exponentials containing x s , which relate to the 

position of the leading edge of the stator along the x axis. For instance, the direct rotor wave 
upstream is analogous to equation 19 above for the stator 

OO OO — N p 

v 1(lir =a£ £ (42) 

n = 1 &=-<*> y'=l 

The direct waves downstream are analogous to equations 20 and 21 and the three coefficients 
corresponding to equation 22 are 
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KR] _ M rl V 1 e -»(«/to)COS0; +3 to sin0 ; )c / z 


o; 




KR2=. 


M 


rl v 2 -«(a 2te c ose ; + P*„sine / )f';Z 0 


( 43 ) 


Pfcn 


05 _ ^_rl_ v 3 ^ -/(a^cose; +p*stae / )c / z 1 
^1 a 3fai 


0/ 


Reflections from the actuator disks are handled with the same reflection coefficients derived for 
the stator waves. Note the permutation of the n's and k's in comparison with equation 26 

VR1 e ia,knXe = (KR2 + VR2 ) e ia2kn X * R 2l (k,n) 

+ (KR3 + VR3 ) e ,a *" x * /? 3 j (k,n) (44) 

VR2 = (KR1 + VR1 )/? 12 (&,n) 

W?J = (KR1 +VRl)R n (k,n) 

where, as for the stator loading case, the VR's can be solved for in terms of the KR’s. 

Effect of Rotor on Rotor 

Upwash caused by the reflected waves on the rotor due to the rotor loading is given by the 
analog of equation 30 for the control points already defined for the stator-on-rotor effect in 
equation 35. 


w 


N„ 


rn 


W, 


= E £ / KRR / (n,iJ)xLR(nJ)e 


inBjili 


(45) 


1 /!= 1 j=\ 

where the array of influence coefficients is given by 

oo 

1 


KRR 1 = Y, 


A— 


M 


rl 


[(p^cosQj -a^sinOj) VR1 e 4 a /to«>se 7 sin e 7 )c 7 z,. 

+ (P^ cos 6 [ - <x 2kn sin 0 j ) VR2 e i< ' a ’ 2kn C0S 6; + si " Q] ^ ]Zi 
+(a 3fc/J cos8 l +p te sin0,) VR3 e l(a ^’ cosQl + Pte s,,10 P^ z . ] 


(46) 


As before, we truncate the infinite Fourier series in equation 43 to the following form in analogy 
with equation 32 


w 


N, 


rn 


W 


= Y WRR'(n,i) e 


inB 2 Qt 


(47) 


1 n=\ 


where 
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( 48 ) 


N n 


WRR ' 0 n,i ) = £ KRR ' ( n,iJ)xLR(nJ ) 

7=1 

When the teims for the effects of the direct wave are added from the original Smith formulas, 
the primes are dropped and equation 46 becomes 


N n 


(49) 


WRR(n,i) = J2KRR(n,ij)xLR(nj) 

7=1 

This is the desired coupling equation for the upwash on the rotor caused by the rotor loading. 


Effect of Rotor on Stator 

Here the waveform, including only a finite number of harmonics, is 


w 


rsi 


N h N p 


w. 


= E E E K R S(n,k,iJ)xLR(nj) e 


ikBjilt 


(50) 


2 n - 1 k=-W A ;=! 


where 


KRS(n,k,iJ) = 

1 .[(^cos0 2 -a^sin Q 2 )VR1 


M 


r2 


KP^cose, -a^sin 0 2 )(*/?2 +VR2)e ia ^° e 


(51) 


- (a 3 ^ cos 0 2 + p^sin 0 2 ) (KR3 +VR3)e 3kn * e 


ia 3k» x s „ i(a Jfa ,cos0 2 +p to ,sine 2 )r 2 z j 


1 


As with the stator loading formulas, we rewrite equation 50 in terms of positive frequencies, 
preserving only the real part of the equation 


w 


rsi 


"h N h N p 


W. 


~ E E E t KRS{n,kJj) x LR{nJ) + [KRS{n, -k,ij) x LR(nj) ] * } e 


m* t J kB I at 


(52) 


2 k=l n= 1 y'=l 

To define WRS(k,i), this can be written 


— = WRS ( k ’0 e lkBl£lt 
W 2 k= 1 


(53) 


where 


N„ Np 

WRS(k,i) = 52 E ( KRS(n,k,ij)xLR(nj) + ( KRS{n, -k,ij) x LR(nJ) ) * } 

n=l 7=1 

which is the last of the coupling equations. 


(54) 
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3.4 SOLVING THE COUPLED EQUATIONS 


Equations 34, 41, 49, and 54 are four sets of linear equations giving the upwash velocities 
at the stator and rotor control points as functions of the unknown loading on the rotor and stator 
panels. The upwash components are considered known from the requirement of flow tangency 
on the blades and vanes. Specifically, the condition to be satisfied on the rotor is 

WR(n,i) = WRR(n,i ) + WSR(n,i ) = 0 , all n,i ( 55 ) 

This states that the sum of the upwash at the rotor due to the rotor and stator waves is zero for 

all harmonics and at all control points. For the stator, the condition is different because it is 
excited by the viscous wake, as indicated in figure 2. If we call the wake upwash harmonics 
WWS(n,i), then the stator boundary condition is 

WS(n,i) = WSS(nJ) + WRS(nJ) = - WWS(nJ), all n,i ( 56 > 

indicating that the net effect of the waves from the stator and rotor must be equal and opposite 
to the wake upwash for flow tangency. Formulas for the wake are given in appendix E. 

The preceding equations can be put into a matrix form written schematically as follows 

KRR KSrtl \lr] [Wk] (57) 

KRS KSSj X [LS\ " [WS_ 

which indicates that the rotor and stator loading each cause upwash at both the rotor and stator. 
To find the loading, the matrix equation is inverted subject to the boundary conditions in 
equations 52 and 53. The method to find the sound pressure and sound power from the unsteady 
loading is described in Section 3.5. 

As indicated above, equation 57 only gives the form of the matrix equation schematically. 
To understand the details, we must look down two levels deeper into the matrix structure. 
Equation 57 shows that rotor and stator are both coupled to themselves and to each other but it 
does not show the frequency scattering explicitly. This can be seen at the next level by 
examining equation 58 where the harmonic dependence is stated explicitly for the special case 
of N h = 3. Loading harmonics are counted by n and upwash harmonics are counted by k. 
From the empty submatrices, it can be seen that there is no frequency scattering for each blade 
row acting on itself. However, for each blade row acting on the other, each loading harmonic 
produces upwash at all harmonics. 

There is a problem which prevents equation 58 from being used directly. Examination of 
equations 4 1 and 54 reveals that two of the coupling equations do not express a simple linear 
relationship from loading to upwash because of the appearance of the loading complex 
conjugates. This is dealt with by expressing the equations in pure real terms and is explained 
below first for the rotor-on-rotor submatrices and then for the rotor-on-stator submatrices. The 
stator loading submatrices then follow immediately. 
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Equation 49 for the rotor-on-rotor coupling can be written in terms of real and imaginary 
parts as follows 

W R +iW i = '£ ( K R + iK l) * ( L R +iL l> 

> (59) 

•J^ K R L R-KiLi)+HK,L R +K R L t ) 

j 

which can be expressed in matrix form 



which is all real, as desired. This breakdown into real and imaginary parts is applied to each 
submatrix of equation 58. If we write the entire coupling system from equation 58 simply as 

J>(p,v) L(v) = W(p) ( 61 ) 

v 


then we place the real parts of KRR into K as follows 

Reai[KRR(n,iJ)] -> K[ (2 n - 2) N +i , (2 n - 2) N +j ] 
and (62) 

Real [KRR(n,iJ)] -> K[ (2 n - 1 ) N p + / , (2 n - 1 ) N +j ] 
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(63) 


and the imaginary parts 

- Imag [KRR(nJJ)] K[ (2 n - 2) N p + / , (2 n - 1 ) N p +j ] 

and 

Imag[/0?R(>w,/)] ^ (2« - 1)A^ +t , (2n -2)^ +; ] 

The corresponding locations for the stator-on-stator elements in the lower right hand comer of 
the K matrix for the real parts are 

Real[KSS(n,ij)] K[ 2 N N h +(2 n -2 )N +i , 2N p N h +(2 n -2)N ' +j ] 
and (64) 

Rca\[KSS(n,ij)] -* K[ 2 N p N h +(2 n - 1 )N p +/ , 2N p N h +(2n -l)N p +j] 

and for the imaginary parts 

- Imag[m(n,iV)l -*K[2N N h +(2n-2)N +i , 2N N h +(2n-l)N +j] 
and (65) 

Imag[m(n,iV)] K[ 2N p N h+ (2n - l)N p +i , 2N p N h +(2n -2)N p *j ] 

The situation for the rotor-on-stator elements is slightly more complicated. First, we rewrite 
equation 54, noting that the conjugate of a product is the product of the conjugates in the second 
term 


"h N p 

WRS(k,i) = £ £ [ KRS(n,k,iJ) x LR(nJ ) + m *(n,~k,ij) x LR *(nj) ] < 66 ) 

n=I y=i 

Then, using the shorthand notation KRS( n, k, ij) =K R (k)+iKj(k), this can be written in terms of real 
and imaginary parts 


N h N p 

WRS{k,i) = £ E [**(*) +*'*/*)] (I* +i'L/) + [**(-*) -i £/-*)](£* -iL,) 

n - 1 7=1 

Separating the real and imaginary parts leads to submatrices of the form 

Ix^ 

\Kfk)~K{,-k)\ [K R (k)-K R (-k)) X 


(67) 


1 


1 

aT 
1 



( 68 ) 


The result of the above development is that the KRS real and imaginary parts are placed into 
the lower left corner of K as follows 


Real [KRS(n,k,iJ)] + Real [KRS(n-k,iJ)] -> K[ 2 N p N h +(2 k -2 )N p +i , (2n -2 )N p +j ] 
-Imag [KRS(n,k,iJ)] ~ Imag [KRS(n-k,ifi] -> K[ 2N p N h *{2k -2 )N p +i, (2 n - l)W p +y ] (69) 
Imag^RS (n,k,ij)] ~ Imag [KRS(n,-k,iJ)] tf[2N p /V A +(2/: - l)N p +i, (2n -2) AT +y] 

Real [KRS(n,k,iJ)] - Real [KRS(n-k,iJ)} *[2^#,, + (2£-l)/V ; , +/, (2n - l)N p +y] 

The development for the KSR terms proceeds exactly the same way to place these elements in 
the upper right corner of the matrix in equation 58. The result is given below 
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Rcal[KSR(n,k,iJ)] + Real[ KSR(n,-k,iJ)] -4 K[(2k -2 )N p +i , 2N p N h +(2n -2)N p +j ] 
-Imag[AT5/?(/i,^, /,/)] - lmag[KSR(n,-kJj)] -> tf[(2*-2)tf p +/ , 2N p N h +(2n -l)tf p +/ ] 
ImagltfS/?^,*,/,/)] - lmag[KSR(n,-k,iJ)} -4 K[(2k - l)N p *i, 2N p N h +{2n -2)N p +j\ 1 ’ 
Real[KSR(nXiJ)] - Real [KSR(n,-k,iJ)] -> K[(2* - l)W p +i, 2W p W A +(2n - l)W p +/] 

Finally, for inversion of equation 58 (or 61), we need to describe the construction of the 
upwash vector W(/j) in terms of the rotor and stator upwash boundary conditions in equations 
55 and 56 and the interpretation of the rotor and stator loading components LR(n,j) and LS(n,j) 
in terms of the loading vector L(v). The upwash vector is constructed as follows 

Real[WK(JU)] -* W[( 2 k- 2 )N p +i] 

Imag [WR(k,i)} -» W[(2k - 1 )N p +i] 

Rea\[WS(k,i)]^W[2N p N h +(2k-2)N p +i] K ’ 

Imag[lV5(it,0] -» W[2N p N h +(2k-l)N p +i] 

and inverting the matrix and computing the L vector, the rotor and stator loads are constructed 
from L(\) according to 

LR(nJ) = L[(2n -2 )N p +j] + i L[{2n - 1 )N p +j] 

LS{nj) = L[2N p N h +(2 n -2 )N p +f) + i L[2N p N h +(2 n - l)lV p +;] 

If we only needed unsteady blade loading, we could stop here. However, for acoustic 
application, we are interested in the upstream and downstream pressure waves and their sound 
power. This is the subject of the next section. 


3.5 SOUND PRESSURE AND SOUND POWER FORMULAS 

The preceding sections showed how to find the unsteady rotor and stator loading in a 
coupled fan stage and gave expressions for the transverse velocity components of the waves 
caused by this loading. For acoustic purposes, we need to derive the expressions for sound 
pressure as well. Furthermore, acoustic power desired as a fundamental measure of sound 
generation independent of time and axial and transverse position in the duct. To develop 
formulas for acoustic power, we also need formulas for the axial component of acoustic velocity. 
The expressions for pressure and axial velocity are developed next and then power formulas are 
developed later in the section. 


Sound Pressure and Axial Velocity 

The theory of this report develops expressions for acoustic disturbances in three regions: 
Region A extends from the rotor leading edge to upstream infinity. Region B lies between the 
rotor leading edge to the stator trailing edge. Region C extends from the stator trailing edge to 
downstream infinity. Section 3 developed expressions for the transverse components of the 
disturbance velocities in Region B. To find the acoustic pressure and velocity in Region A, the 
required steps are to 1) find the transverse velocity at the upstream boundary of Region B, 2) 
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apply the transmission coefficients developed in appendix D to jump to the transverse velocity 
in Region A, 3) apply Smith’s equations (equations 6 and 7 herein) relating pressure and axial 
velocity to transverse velocity to find the desired quantities. Acoustic pressure and velocity in 
Region C are found by the same method. 

Waves in Upstream Region 

We start with the stator waves given by equation 27. We select only the upstream-going 
waves since definition of the transmission coefficients is based on them. At the rotor leading 


edge, x=0 + , 


oo oo 


AL 


(73) 


vj « . a b £ £ j * iyVnk P nk (KSl e + VS1 ) x LS(nJ) 

71 = 1 fc=-°o y' = l 

To find the corresponding value just across the inlet actuator disk at J=0 , we apply the 
definition of the transmission coefficient in equation 16 (and in appendix D) to find 


oo oo 


N n 


v 4 =a aE E E « Kk( KS1 e ~ iainkX ° + VS1 ) T \aM) * ls m 

71=1 k= )=1 

Now we define 

KSUP(n,kJ) = (KS7 6! + VS7 ) r 14 (n,*) 

and insert the axial wavenumber exponential to find v 4 at any upstream location 


oo oo 


AL 


h=1 k~ y=l 

The working form for this is obtained by performing the y'-sum and defining 


(74) 


(75) 


(76) 




LSUP{n,k) = 53 KSUP(n,kj) x LS(»7) 
7=1 


(77) 


so that 


oo oo 


^^E E 


, «(V»* + a***) 


LSUP(n,k) 


(78) 


71 = 1 fe=-oo 


As mentioned above, equations 6 and 7 provide the connection from the transverse velocity 
component to the axial component and the pressure. When these are applied, we find 


OO oo 


stator 


= a aH E <*. 

fl=i k - -« 


4/2 k 


, ‘(Vnk + a 4nk x ) 


LSUP(n,k) 


(79) 


and 
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OO CO 


P 4 


stator 


= ?oa a l E E 


where we have defined 


n=I <:=-<» 


fli 


Tlrcfc ( + ^ xa^Ank) 

a„ 


(80) 


( 81 ) 


Equations 79 and 80 are the desired forms for the axial velocity and pressure in the upstream 
region caused by stator loading. Derivation of the corresponding forms for waves from rotor 
loading is similar and leads to the following results. 


and 


oo oo 


i rotor 


■«,EE a 4ta / (¥ ‘"‘‘ W) LRUP(n.k) 


n-l k~-°° 


(82) 


\Pa 


rotor 


m Poa a a E E y ikn e‘^' a ^ X) LRUP(,n,k) 


n = 1 k=- 


(83) 


Permutation of the n and k indices (in comparison with equations 79 and 80) arises from 
viewing the rotor waves in the stator coordinate system. This is an important aspect of the 
kinematics of blade row interaction and is discussed in detail in appendix A. The definitions of 
LRUP and KRUP are similar to those for LSUP and KSUP above: 


and 


LRUP(n,k) = 52 
7=1 


KRUP(n,kJ) x LR(nJ) 


KRUP(n,kj) = (KR1 +VR1 ) T l4 (k,n ) 


(84) 


(85) 


The total acoustic pressure in the upstream region, of course, is the sum of equations 80 and 83. 
By noting that appears in equation 75 but not in equation 85, it can be seen that the phase 
difference between the rotor and stator waves associated with axial spacing is accounted for and 
thus that interference effects associated with rotor/stator spacing are included in the theory. 


Waves in Downstream Region 

This derivation is similar to the one just presented, but slightly more complicated because 
the jump between regions is not at x=0 and because both pressure and vorticity waves impinge 
on the downstream actuator disk. From equation 28, the downstream-going waves from the stator 
are 


oo oo 


AL 


;*v,=a„ E E E 

n-[ k~ j= 1 _ _ 

+ a 3nk (KS3e + VS3) e ] x LS(nJ) 


‘ a 2nk x s + YS2) P 


( 86 ) 
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at the exit plane x = x e '. The definitions for the transmissions coefficients T 28 and T 38 from 
Equations 17 and 18 are applied to find the transverse velocity in the pressure wave at x = T e + 
and the x-wavenumber exponentials are added to find the velocity anywhere downstream of the 
stator actuator disk. 


E \ k W2 e -‘ a ™*- + VS2)e ia ^‘T n (n,k) 

n- 1 k- -<*> 7=1 

+ (KS3 e ~ ia ***' + VS3) e r 38 (njfc)] e ia8nk(x ~ Xe) xLS(nJ) 

Now, by defining 

KSDN(n,kJ ) = (KS2e + KS2)/ a2n ^ ^(n,*) 

+ (KS5 e + VS3) ^ ' “***« r 38 (n,it)] 

and 


(87) 


( 88 ) 


n p 

LSDN(n,k ) = £ KSDN(n,kJ) x LS(« j) 

7=1 

the transverse velocity for the pressure wave in the downstream region can be written 


v 8=«cE E Pnfc « ' IV "* +asni:(j: ~ Xe)] LSDN(n,k) (90) 

n= 1 k- 

Again by applying equations 7 and 8, the axial component of the acoustic velocity and the 
acoustic pressure in the downstream region are found to be 


and 


oo oo 


1 stator 


= a c E E a : 


8 nk 




LSDN{n,k ) 


n-i k~ -<» 


P 8 




stator 


Por a c E E "^8 nk e 


LSDN(n,k ) 


«=1 k - -oo 


where y gnk is given by 


nk ~ 


^nk + ^xc Ck Snk 


( 91 ) 


(92) 


(93) 


For the rotor waves the procedure is similar and leads to the following forms for axial 
velocity and pressure 


22 



(94) 


and 


©o oo 


[“«Loc* a ‘£ E ^ LRDNM 


n = 1 k=-oo 


where 


oo oo 


hLr'^E E 


n-\ 


(95) 


and 


LRDN(n,k ) = £ KRDN(n,kJ) x LflOy) 
7=1 


KRDN{n,kJ) = (KR2 +VS2)e‘ a2k * X ° T 28 (k,n ) 

+ (03 + VR3) e ' T 38 (j(,rt) 


(96) 


(97) 


The above equations for pressure are coded in CUP2D for output on a modal basis and as 
a sum over the modes at each frequency. The pressure formulas and those for the axial velocity 
components are used in the next section to derive equations for acoustic power. 


Sound Power 

Formulas for sound power are derived from Goldstein’s (ref. 9) equation for acoustic energy 
flux vector 

7=(p/p 0 +u-v 0 )(p ( ,«+pv 0 ) (98) 

where p , p , and u are the acoustic pressure, density, and velocity and v 0 is the mean 
velocity of the background flow. We are interested in Regions A and B, where the flow is axial, 
so that u v 0 =uU. Then, with p -p/a 2 , the flux becomes 

/=(p/ P<> *t/ U )(p„<, (99) 

a 1 

Since we are interested in the flux in the axial direction, we form 

l x =Ii x = ( plp 0 + Uu) (P 0 U + JL U) (100) 

a 2 


or 
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(101) 


I x = (1 +M 2 x )pu + —Lp 2 + p 0 aM x u 2 

This is the instantaneous energy per unit time crossing a unit area element that is normal to the 
x-axis. The power is the time average of I x integrated over the duct cross-sectional area. Since 
the concept of area in a 2D duct is nebulous, we will simply compute a non-dimensional quantity 
equal to the average power per unit area normalized by the ambient value of p 0 a 3 . The average 
power is then 

(102) 

Since this is defined to be positive for flux in the x direction, computed values for our problem 
must be positive downstream of the source and negative upstream to represent energy propagating 
away from the source region. 

Upstream Power 

To begin operations with equations 101 and 102, we need convenient forms for u and 
p. As stated above, the axial velocity in Region A is given by the sum of the stator and rotor 
contributions expressed by equations 79 and 82. 

u A = a a Real£ £ a Anlc LSUP(n,k) +a4nkX) 

"= 1 ~ (103) 

+a a Real E a 4ta LRUP(n,k) e i(v *" +w) 
n = 1 fc=-«o 

The declaration "Real" has been implicit all along but is shown here explicitly because of the 
need to eliminate it for the manipulations that follow. This is done by changing to double sided 
series in n as follows 


1 n=- °° k=-°° (104) 

+ — E E a Akn LRUP(n,k)e i(Whl+a4knX) 

2 n - -oo k=-°° 

The reader can verify that equations 103 and 104 are equivalent by noting that and o. 4nji 
simply change sign when both indices change sign and by defining LSUP(-n,-k)=LSUP (n,k). 
The "Real" declaration in no longer needed since the imaginary parts of the upper and lower 
halves of the n-series cancel identically. The key step now is to switch the dummy k and n 
indices in the second line of the equation, yielding 

« 4 = — E E a 4nk [LSUP(n,k) + LRUP(k y n)] e i(Vn * + ( 105 > 

2 n -~ oo k=-°° 

or, by defining 
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LUP(n,k ) = LSUP(n,k ) + LRUP(k,n) ( ■ 106 ) 

we obtain the simpler form 

“4 = -^- E E ( 107 > 

^ n =-° o - oo 

which will be seen to be convenient for manipulation later. 


By a similar argument, the desired form for the upstream pressure is 

E y^LUPinJOe^'*^ 

^ /J=-oo /c= -OO 


(108) 


Now, consider the expression for the product p 4 u 4 =p 4 u 4 (because is real) which has be 
averaged over (|) and t. This requires the quadruple summation 


Pa u a 


.1 

—7— EEEE ««T4 mq LUP(n,k)LUP\m,q)e i{a4 " t ~ a - 4 ^ )I 


n m k q 


(109) 


-i[(n -m)B] -(k.-q)B 2 ]ty i[(n -m)B 1 Cl l +(k -q)\B 2 Q. 2 \t 

X ^ (f 


Averaging over (() involves only the second exponential. That average is zero for each term in 
the series except for the cases where (n-m)Bj = (k-q)B 2 , when the average value is unity. When 
this expression is substituted into the time exponential in equation 109, the result is exp[/(/z- 
m)B 1 (Q. 1 +Q. 2 )t]. When this is averaged over time, the result is unity for n=m and zero 
otherwise. Thus, when the m and q sums are performed in equation 109, the exponentials 
vanish and the result is 


T, 

^ C° CO 

— = £ a tnk y inll LUP(n,k)LUP\nJc) < 110 > 

^ n = — oo oo 

where the double overline denotes the average over <t> and t. From the symmetry properties 
of the a’s, y’s, and LUP' s, it can be shown upper and lower halves of the /i-series are identical 
so that 


— £ <W4 nk LUP(.n,k)LUP'(nje> OH) 

^ n= 1 Jt=-o° 

where the n=0 term is omitted because it makes no noise. This is the working form for the 
average of the first term in equation 101. It exhibits the usual behavior that the total power is 
the sum of the modal powers. Also, as required, the result is independent of 3c and <|). 
Averages of the remaining terms in equation 101 can now be written down by inspection. The 
final form for the sound power in Region A is 

E I(l+<)a«T«^„(a^ + yL)]l i -WW)l 2 (ll2) 

^ n-\ &=-oo 


25 



This is the form used in the computer code except that the factor '/ 2 p oa a a 3 is omitted via the 
following normalization. 


n xa = 


n 


x a 


■tP oa a a 


(113) 


For a decibel output, the value is divided by 10' 13 , arbitrarily, to give positive values for plotting 
purposes 


PWL = 10 log 10 


n 

10 


xa 


(114) 


Downstream Power 

The derivation for power in Region C follows along identical lines. The result is 


P oc a c 


n xc = E K 1 + ^)«8 nk fs nk +M x c«hnk + ink^ I U)N{nJC) | 2 (U5) 

^ n = 1 k.=-°° 


where 

LDN(n,k) = LSDN(n,k) + LRDN(k,n ) (ll6) 

Normalization for this in the code is the same as for the upstream power, except that the 
downstream conditions and a c are used. 


The theory of this section is coded in program CUP2D. Documentation for CUP2D can 
be found in Volume 2 of this report and includes a source listing, user instructions, and a test 
case. Typical cases for 3 harmonics and 30 panels on each blade row run a few minutes on a 
workstation computer or a 386 PC. 
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SECTION 4 

COMPARISON WITH OTHER THEORIES 


The theory of the present paper adds three effects not included in the Smith theory: blade 
row coupling, frequency scattering, and flow turning. For verification of the theory, it was 
possible to check only two of these effects by comparison with previous work. 

Comparison with Kaji and Okazaki 

Kaji and Okazaki developed a method (ref. 6) to compute the coupled interaction between 
a downstream rotor and an inlet guide vane (IGV) cascade with viscous wakes. Their method 
could treat one harmonic only and did not account for turning of the flow. Since the present 
analysis treats rotor/EGV, rotor/rotor, and IGV/rotor interaction, it was possible to make 
comparisons with the earlier Kaji/Okazaki results. (The wake model used here is equivalent to 
that of Kaji and Okazaki and is developed in appendix E.) A sample is shown in figure 8 for 
a case where B 1 -B 2 and the gap/chord ratio is the same for both rotor and stator. The 
geometry and flow condition are such that the relative Mach number is the same for both blade 
rows. The curves show the upstream and downstream sound pressure for two rotor spacings as 
relative Mach number varies. Agreement of the two theories is excellent. This provides 
confidence in the coupling part of the new theory, which is handled differently in the Kaji and 
Okazaki work. 

Simulation of the Kaji/Okazaki case was as follows. The IGV and rotor blade angles were 
both 30 deg., gap/chord ratios were both 1.0, and the interblade phase angle was zero. Since 
theirs was an unloaded model, this corresponds to a stage with an inlet (and exit) flow angle of 
30 deg. The model of this report does not permit swirl in the exit flow and therefore the Kaji 
and Okazaki configuration was simulated by running both blade rows as rotors with equal and 
opposite RPMs. By using -M y2 = M y] = M x tan30 deg. , the geometry and flow angles were 
correctly represented. The blade numbers had to be equal and were chosen to be equal to 10. 
Conditions were taken to be sea level standard for density and sound speed. 

Comparison with Kousen and Verdon 

For single cascades, a much more sophisticated method has been developed at the United 
Technologies Research Center by Verdon and his co-workers (ref. 10). Their code, called 
LINFLO, starts with a steady flow field from a compressible potential flow calculation and 
computes an unsteady perturbation field due to input waves or blade motions based on a 
linearized, inviscid, unsteady analysis. In a recent paper, Kousen and Verdon (ref. 1 1) showed 
calculations that provided an opportunity to establish how well the present method represents the 
effect of flow turning. In figure 9 the cascade is excited by a vorticity wave with an interblade 
phase angle and frequency such that it is cut on upstream of the cascade and cut off downstream. 
Kousen and Verdon compared this to a no turning, flat plate calculation (equivalent to a Smith 
code calculation) and found that the case with turning increased the upstream response wave by 
a factor of 8.4, or 18.5 dB. Some alterations to the present method were made to represent a 
single flat plate cascade with (or without) turning at the exit. When the Kousen/Verdon case was 
simulated, it was found that the reflected wave was amplified by 19.7 dB by the turning. It was 
concluded for this case that the primary effect of the real blades (beyond the flat plate results) 
is reflection of the waves from a region where they cannot propagate and that the simplistic 
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actuator disk model did a reasonably good job of representing the effect. 

Simulation of the Kousen/Verdon case was as follows. In the upstream/downstream regions 
the axial Mach numbers were 0.476/0.450, the swirl Mach numbers were 0.275/0.000, the sound 
speeds were 1 1 16/1 128 ft/sec, and the densities were 0.07631/0.080185 lbm/ft 3 . Use of B l =36, 
B 2 =57, chord/radius = 0.131, and rotor rotational Mach number = 0.478 produced the correct 
reduced frequency (4.106) and interblade phase angle (2.31). To eliminate wave reflections from 
the rotor and inlet actuator disk, KR1, KR2, KR3, R J2 , and R J3 were set to zero. Finally, 
reflections from the exit actuator disk (representing the mean flow turning) could be turned off 
by setting R 21 and R 31 equal to zero. 
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SECTION 5 

BEHAVIOR OF THE COUPLED CASCADE THEORY 
& THE MODE TRAPPING PHENOMENON 


In this section, behavior of the new theory is examined using a format similar to that of the 
test data in figure 3 to find to what extent the features of coupling, frequency scattering, and flow 
turning explain the observed behavior. Since the theory is 2-dimensional, we are only looking 
for qualitative behavior at this point; absolute level predictions will have to await a 3D analysis. 
The calculations lead to the conclusion that the high levels of noise at 2x and 3xBPF are caused 
by a mode that reflects back and forth between the rotor and stator, amplifying in level and 
coupling to the upper harmonics. The final part of this section verifies the assertion made during 
derivation of the theory that valid results can be obtained by including only a small number of 
loading harmonics in the coupling equations. 

The cases presented in this section are all in the form of RPM sweeps with mean flow 
conditions representative of the 85 percent radius in the turbofan of figure 3. The mean flow 
values used for the computations, which are given in Table I, were obtained by running a 
standard streamline curvature code. Each column corresponds to one operating condition. For 
the curves to be presented below, M yJ (the first row) is treated as the independent variable and 
values of the other parameters are found by straight line interpolation on M y] . 


TABLE I - MEAN FLOW CONDITIONS 
a) Upstream of rotor, b) Between rotor and stator, and c) Downstream of stator 


im 

0.672 

0.908 

0.974 

1.017 

M xa 

0.269 

0.418 

0.467 

0.496 

M xb 

0.321 

0.440 

0.469 

0.480 

M xc 

0.382 

0.519 

0.544 

0.546 


0.205 

0.285 

0.334 

0.359 

P oa 

.02978 

.02837 

.02767 

.02738 

Pob 

.03222 

.03407 

.03525 

.03592 

P oc 

.03201 

.03381 

.03539 

.03647 

a a 

1041 

| 

1031 

1027 

1025 

a b 

1065 

1080 

1091 

1098 

a c. 

1065 

1081 

1094 

1104 


p 0 - density - lbm/ft 3 , a - speed of sound - ft/sec 
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In particular, the reflection and transmission coefficients are computed from the jumps in mean 
flow properties from regions a to b (across the rotor) and from regions b to c (across the stator). 
For comparison with Figure 3, note that M yl is taken as 85 percent of the fan tip rotational 
Mach number. In addition to conditions given by the table, the rotor and stator gap/chord ratios 
were 0.763 and 0.957 and the spacing from the rotor trailing edge to the stator leading edge was 
two rotor chords. There were 38 blades and 72 vanes. Stagger angles follow the mean flow 
according to Figure 5 but were generally about 55 deg. for the rotor and 32 deg. for the stator. 
Furthermore, all were run with a wake model equivalent to that in the Kaji/Okazaki work 
(derived herein in appendix E) with a drag coefFicient of 0.02 but with all of the wake harmonics 
above BPF set to 0. Thus, excitation of the coupled system is exclusively at BPF and any 
appearance of higher harmonics must be caused by frequency scattering. These conditions apply 
to all cases discussed in this section. 


5.1 EFFECTS OF COUPLING, FREQUENCY SCATTERING, AND SWIRL 

Four RPM sweeps are shown in figure 10 to demonstrate the effects of the new theoretical 
features. The plot at the upper left presents a baseline case and then various effects are "turned 
on" to illustrate their influence. For the cases without swirl, all the mean flow conditions are set 
to the values for the aft duct in the table above. Thus, without swirl, all of the reflection 
coefficients associated with turning of the mean flow are 0. 

Figure 10a shows calculations without swirl for one harmonic only. BPF cuton in the aft 
duct is at a rotor rotational Mach number of 0.80. The first mode to cut on at BPF has order 
m=- 34, which is computed from m=nB 1 -kB 2 with B ; =38, B 2 =12 , n- 1, k= 1. The uncoupled 
curve corresponds to results from the original Smith theory (KRS and KSR were set to zero). 
The case with coupling uses theory equivalent to that of Kaji and Okazaki. It can be seen that 
the coupling effect is very weak for this case and could easily be ignored for acoustic purposes. 

Figure 10b was generated without swirl but with coupling and scattering to 2 and 3 times 
BPF included. This is the first case showing excitation at BPF producing noise at higher 
harmonics. Also, note that there is a small amount of noise in the higher harmonics at speeds 
below BPF cuton. This figure does not look at all like figure 3 with its high levels of 2xBPF 
and 3xBPF below BPF cuton. However, above cuton, the scattering effect is not negligible. For 
example, a spectrum at a rotor rotational Mach number M y] = 0.82 would show 2xBPF and 
3xBPF only 5 dB below BPF. 

Figure 10c shows coupled calculations for one harmonic only, with and without swirl. The 
change here is due mostly to the fact that, in the theoretical model, the cascade stagger angles 
follow the mean flow as indicated by figure 5. Thus, for the case without swirl, the stator angle 
is 0 deg. and the rotor stagger is typically 60 deg. The extended range of rotational speeds with 
swirl is due to the fact that swirl reduces the relative Mach number at the rotor fig. 5) and 
permits a higher RPM before the rotor relative Mach number reaches unity. The rapid rise in 
the curve for the "with swirl" results from cuton of a second propagating mode ( n,k = 1,0) at 
M y] =0.87. 
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Finally, figure lOd shows calculations with all of the theoretical features included. We see 
that the combination of swirl and coupling produces a substantial increase in levels of the higher 
harmonics, particularly below BPF cuton. In comparing this with the plot of engine data in 
figure 3, we see that, although the plots are not exactly the same, there are some striking 
similarities. At the RPM where BPF cuts on in the interstage area, 2xBPF and 3xBPF rise 
rapidly and are actually higher than BPF when it cuts on in the aft duct. When BPF cuts on in 
the aft duct, the higher harmonics drop substantially. In the intermediate speed range, 3xBPF 
is higher than 2xBPF by about the same amount as in the engine data. These similarities suggest 
that the 2D model embodies the correct physical features to explain the experimental behavior. 


5.2 INTERPRETATION OF RESULTS 

To interpret the results presented in figure 10, it is necessary to understand behavior of the 
reflection coefficients associated with the inlet and exit actuator disks and how they are related 
to modal cutoff ratios. The 1,1 mode is the first to cut on and we show its behavior in figure 
11. The upper panel relates to the inlet interface and shows R I2 (I,1) and R ]3 (],J), which are 
the coefficients for an upstream going pressure wave reflecting into downstream going pressure 
and vorticity waves. Point b indicates cuton between rotor and stator and point a indicates 
cuton in the inlet duct. It can be seen that reflection into the vorticity wave is always weak. 
However, for the pressure wave, the coefficient is of order unity between a and b and is small 
outside. This indicates that a wave travelling from a region where it is cut on to a region where 
it is cut off is reflected with little change in amplitude. The lower panel in figure 1 1 relates to 
the exit interface and shows R 2J (1,1) and R 3I (1,1), the coefficients for the upstream going 
pressure wave reflected from the downstream going pressure and vorticity waves. Again, there 
is strong reflection for waves traveling from a region where they are cut on to one where they 
are cut off. This is believed to be the dominant reason for the strong reflections shown in the 
Kousen/Verdon^ 11 ^ calculations in figure 9. 

The explanation for the behavior of the upper harmonics shown in figure lOd can be 
explained pictorially with the sketch in figure 12. The top panel relates to the speed range from 
cuton between blade rows to cuton in the aft duct. In this range, the 1,1 mode can bounce back 
and forth freely between rotor and stator, but it cannot escape. It effectively is trapped in the 
interstage area and builds to high levels as it couples to the higher harmonics. The higher 
harmonics then can escape to the aft duct. Thus, in this speed range, the excitation of the system 
is BPF of the rotor wakes, but the only noise escaping from the system is at higher harmonics. 
The lower panel in figure 12 relates to speeds where BPF can propagate in the aft duct. The 
1,1 mode is no longer trapped so that the BPF levels between rotor and stator do not amplify 
as strongly and the coupling to the higher harmonics produces lower levels. 

Finally, note that figure 10 shows that the combination of coupling and swirl was required 
to produce the effect. Swirl or coupling by itself produced much less dramatic results. 
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5.3 CONVERGENCE WITH INCREASING NUMBER OF HARMONICS 


In the theory derivation of Section 3, the coupling equations were initially written as waveforms 
of upwash with an infinite number of harmonics. Since the coupled system of equations is 
solved numerically, How tangency can be guaranteed simultaneously at only a finite number of 
harmonics. Thus, the summation over loading harmonic order was truncated at N h rather than 
running to infinity. Since each loading harmonic on the stator produces an infinite number of 
upwash harmonics on the rotor and vice versa, it was not known a priori that this truncation 
would produce acceptable results. The only way to verify correct behavior was to run a series 
of calculations with increasing values of N h . It was hoped that as more and more harmonics are 
included in the coupling, the influence on the lower harmonics would be minimal. 

Verification of this convergence is shown in figure 13, which is based on the same 
conditions used for figure 10. The top panel shows predicted aft PWL at BPF for separate 
calculations where 1, 2, 3, 4, and 5 harmonics were included in the coupling. It can be seen that 
near cuton the number of harmonics does influence results but that levels tend to converge for 
larger N h . When the next mode cuts on near rotational Mach number M yl = 0.87, the 
sensitivity to number of harmonics disappears. It seems that the upper harmonics are taking 
energy from only the m = 38-72 = -34 mode. 

The second panel from the top in figure 13 shows results for 2xBPF. Here the sensitivity 
to N ) j is minimal. The lower panels in the figure indicate that the situation improves at higher 
harmonics. Thus, we have given an example where results for the lower harmonics are not 
strongly affected by the number of harmonics included in the coupling. Although all of the cases 
checked so far exhibit this behavior, this does not constitute proof of convergence. A probable 
explanation for the convergence is that loading response diminishes with increasing frequency. 
Also, interference between waves from adjacent panels increases with increasing frequency. 
However, until a more rigorous theoretical explanation can be provided, we cannot be sure that 
convergence will always occur. 
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SECTION 6 

PARAMETRIC STUDIES 


The calculations of the previous section were chosen to provide a physical understanding 
of noise generation in a coupled environment. Geometry and operating conditions were the same 
for all of the figures and are representative of modern turbofan design; in particular, there were 
more vanes than blades. In this section, we explore the new theory via additional computations. 

First, in figure 14 we compare upstream and downstream harmonic powers from the same 
calculation. The top panel in the figure is the same as in figure lOd and is the "base case" for 
variations in the remaining figures. The forward PWL’s shown in the lower panel exhibit 
behavior similar to the aft noise but with cuton of BPF at a slightly higher rotor speed because 
of differences in inlet and exit mean flow conditions. Note that the forward noise levels are 
much lower than the aft noise. This is believed to be due to the fact that the basic noise source 
is the stator and its noise must propagate upstream through the rotor. Note also that the forward 
levels are so low that they are likely to be masked by other sources, such as inlet distortion/rotor 
interaction in a real engine environment. 

The top part of figure 14 also indicates the cuton speed for relevant modes via the arrows 
along the bottom grid line. For example, the notation 1,1b at rotational Mach number M } 
= 0.665 indicates that the n,k =1,1 mode cuts on at this speed in region b (between rotor and 
stator). In the figure it can also be seen that the same mode doesn’t cut on in the aft duct (region 
c) until My] reaches 0.800. Thus, the 1,1 mode is trapped between these rotor speeds. The 
figure also shows that mode 1,0 cuts on at M y , = 0.870, accounting for the rapid rise in BPF 
at this speed. 

Because of the motivation to understand the importance of frequency scattering, all of the 
computed noise so far in this report has been excited by the BPF harmonic (fundamental) of the 
rotor wake; the upper wake harmonics have been set to zero. In the lower portion of figure 15, 
the wake harmonics are computed in the normal way from the wake formulas in appendix E. 
Thus, there are three wake harmonics and three harmonics in the coupling calculation. By 
comparing with the base case in the upper portion of the figure, it can be seen that BPF and 
3xBPF are only slightly affected, however, 2xBPF changes considerably. Since this is a linear 
analysis, superposition applies: the 2xBPF oscillations must be caused by interference between 
noise caused by the wake upper harmonics and noise caused by scattering from the fundamental. 
Note that the overall noise is not increased much by including the upper wake harmonics; that 
is, the upper noise harmonics are caused primarily by the scattering phenomenon and not by the 
better known process whereby each wake harmonic couples only to the same noise harmonic. 
This behavior, of course, depends on the wake properties so that no universal conclusions can 
be drawn. 

The cases above were all for geometry with 38 rotor blades and 72 stator vanes, that is for 
a typical cut off design. For these, the first mode to cut on at BPF has order m=-34, (computed 
from m=nB I -kB :> with /j=l, for BPF, and k=l, for the first wake spatial harmonic). The 
significance of the negative mode order is that this mode spins counter to the rotor direction, and 
the swirl direction, and therefore cuts on earlier in the swirl region than in the aft region. Thus, 
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the typical behavior that we have seen above is high levels of 2xBPF and 3xBPF when BPF is 
cut off. In figure 16, we compare this base case for 38 blades and 72 vanes with calculations 
for the same blade count but 38 and 26 vanes. Blade row solidity and rotor/stator spacing are 
the same for all three vane counts. For fi 2 =38 and 26, BPF is cut on throughout the speed range 
of interest and we do not find any speed range where the upper harmonics dominate, however, 
there is considerable frequency scattering. (All portions of figure 16 were computed with only 
the wake fundamental harmonic for excitation.) Now, for the 38-vane case, note the rapid rise 
in the 3xBPF harmonic near M y] =0.12. This is caused by the 1,2 mode which cuts on at a 
rotor rotational Mach number of 0.72 between blade rows and at 0.87 in the aft region. As 
denoted by the *’s above the 1,2 b and 1,2 c notations, this is another trapped mode. Here 
is a case where BPF noise is cut on, for k= 1, but there is considerable influence on a higher 
harmonic due to another BPF mode which is cut off and counter rotating. For the 26-vane case 
at the bottom of figure 16, similar remarks apply for the 2xBPF and 3xBPF harmonics. In this 
case the 1,3 mode is trapped between rotor speeds 0.750 and 0.915. Thus, we find that 
frequency scattering and mode trapping can have a significant impact on the noise even for cut 
on designs. 

The final set of calculations is for variations in rotor/stator spacing. Figure 17 compares 
aft power levels for axial spacings of 1.0, 1.5, and 2.0 rotor chords between the rotor trailing 
edge and the stator leading edge. There is very little difference in the three sets of curves except 
for a slight decrease for increased spacing. Since all of the curves go down as a family, this 
appears to be due simply to the weaker wake excitation as spacing increases. (The wake 
formulas are given in appendix E.) At the low speed end of the curves, note that the slope is 
steeper for the greater spacing. Noise at rotor speeds below 0.675 rotational Mach number is 
caused by excitation of the rotor by a cut off mode from the stator. Thus, at a given RPM in this 
low range, noise will depend strongly on rotor/stator spacing because of the exponential decay 
of the cut off mode. 


34 


SECTION 7 

CONCLUDING REMARKS 


An analytical model has been developed for rotor/stator interaction by extending the 
industry-standard 2-dimensional unsteady cascade theory of S. N. Smith to include 1) unsteady 
coupling of rotor and stator, 2) frequency scattering due to relative motion between blade rows, 
and 3) turning of the flow due to the mean loading on the blade rows. Unsteady loading is 
handled rigorously by collocation at many panels simultaneously on each cascade for several 
harmonics. The effect of finite mean loading is treated simplistically via actuator discs that turn 
the flow at the rotor leading edge and straighten it at the stator trailing edge. For typical 
turbofans designed for cutoff, there is a speed range where the fundamental mode at BPF can 
bounce back and forth between the rotor and stator but cannot escape to the inlet or exit ducts. 
In this range, the mode amplifies and couples to the higher harmonics, producing a significant 
amount of noise at these frequencies. This mode trapping mechanism appears to explain 
observations in engine data where the source of excitation is at blade passing frequency but the 
spectrum is dominated by 2xBPF and 3xBPF. 

The computational method requires truncating the Fourier series that represents the blade 
and vane loading. It was demonstrated that, when three or four harmonics are included in the 
coupling, noise at the lower harmonics is predicted with reasonable accuracy. 

Limited verification of the theory was possible by comparison with earlier work of Kaji and 
Okazaki, although their analysis did not include frequency scattering or turning of the flow. By 
comparison with calculations from the LINFLO code for loaded, isolated cascades, it appears 
that turning of the mean flow is the major feature of real cascades (as opposed to flat plates) that 
influences noise. 

Although unsteady coupling produces frequency scattering under any conditions, swirl is 
the critical ingredient that produces the mode trapping and amplification. Scattering into the 
upper harmonics can be significant for both cut off and cut on designs. It is concluded that any 
further modeling of rotor/stator interaction for noise calculations should include the features of 
coupling, frequency scattering, and flow turning. 
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APPENDIX A 

KINEMATICS OF ROTOR/STATOR INTERACTION 


The objectives of this appendix are 1) to find, via a kinematic analysis, the form of all of the waves 
that can result from rotor/stator interaction, 2) to relate the interblade phase angle notation of cascade 
theory to the spinning mode order notation of duct acoustics, and 3) to develop one system of indices for 
the transverse wavenumbers that can be used for both rotor and stator waves. 

Kinematics of the interaction are determined by a part of the exponential in equation 1 which we call 
the kinematic phase 

V|/=py+OH 


In the following, we study the properties of y for the rotor and stator and relate it to the spinning mode 
notation of fan acoustics. 

Consider the two blade rows shown in figure 5. For generality, it is assumed that both blade rows can 
rotate in either direction. To simplify the nomenclature, the front blade row will be called the "rotor" and 
the rear blade row will be called the "stator"; however, in reality, the analysis can treat IGV/rotor and 
rotor/rotor configurations as well as the rotor/EGV configuration indicated by the nomenclature. There 
are three coordinate systems with common values for x : the x,y system is fixed in the duct; the x,y 1 
system is locked to the rotor; and the x,y 2 system is locked to the stator. The rotor is shown with 
positive stagger angle 0 2 in conformity with Smith’s nomenclature, but for the rotor/EGV geometry 
emphasized in this report, the rotor stagger angle will be negative. Blades are aligned with the mean flow 
direction so that stagger angles are determined by the axial velocity U the swirl velocity V s and the 
angular speeds according to the velocity triangles in figure 5. The constant R used below can be 
considered the effective radius of a thin annular duct and is the reference used throughout the paper for 
non-dimensionalizing wavenumbers, chords, and distances. The y’s and the cylindrical coordinate angles 
are related by 

y = R<\> 

yi=/?<t>i (A-2) 

y 2 =/?<t> 2 

The coordinate systems are related to each other by 

(|)=(|)2+Q2 ? (A-3) 

<J > 2 = <)>i t 

where 

Q = Qj -ft 2 (A-4) 

Because the Smith equations were developed in a coordinate system locked to a blade row, the 
discussion below will jump from stator coordinates to rotor coordinates and back again. Implications for 
frequencies in the duct-fixed system will be discussed afterwards. Starting in the stator frame, we can 
see that the relative motion of the rotor with B l blades will induce loading harmonics of order nB l with 
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(A-5) 


frequencies nBjQ. The resulting kinematic phase can thus be written 

Inserting equation 2, the relation between blade gap and blade number, S 2 =2kR/B 2 , and the conversion 
to cylindrical coordinates via equation A-2 leads to 

V| f s = (—B 2 + kB 2 )(\ >2 +nB l Qt (A-6) 

2k 

Now, to find the interblade phase angle for the stator o s , compare the condition just analyzed with that 
when the rotor has advanced one stator vane gap, that is by angle 2 kJB 2 in time 2k/B 2 Q. It can be seen 
that the geometry will be identical to figure 5 but with the y 2 (or (}>,) axis shifted. We conclude that \\f $ , 
with <J> 2 — > <J) 2 2 nlB 2 and t — > t + 2 kIB 2 Q., is the same as \|/ v in equation A-6. By equating these 

expressions for \|/ 5 , we find 

a s = -2KnB l /B 2 (A-7) 

for the interblade phase angle, that is the difference in loading phase on two adjacent stator vanes at 
harmonic n. Inserting this into equation A-6 gives the expression familiar from duct acoustics 

= -{nB j -kB 2 )§ 2 +n Bfiti (A-8) 

which indicates modes of order nBj-kB 2 (for all k) appearing at frequency nBjQ and spinning at speed 
[nB 1 /(nB 1 -kB 2 )]Q in stator coordinates. 

A similar derivation for rotor loading harmonics of order n and frequency nB 2 Q leads to interblade 
phase angle 


o r = +2nnB 2 /B l (A-9) 

where the + sign appears because the relative motion of the stator is opposite that of rotor. With these 
results, the kinematic phases for rotor loading harmonic n are 


\\f r = -(kB j -nB 2 )<t>i + nB 2 D.t 


(A- 10) 


for all k. 

We now define the mode order notation for the stator 
and take advantage of the symmetry with the rotor waves to note that 


(A- 11) 


(A- 12) 
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(A-13) 


Then the kinematic phases for the stator (in the stator coordinate system) can be written 

2 +nB \ Qt 

and for the rotor (in the rotor coordinate system) 

\\f r = (3^(1), +nB 2 Qt (A- 14) 


A revealing step now is to find the rotor wave properties in the stator coordinates by substituting <J)^ 
= <f>, - Q.t into equation A- 14. The result is 

Vr = Mi +kB i Qt (A- 15) 


which is to say 


Vr = V*« 


(A- 16) 


Since indices n and k take on all integer values, comparison of equations A-13 and A- 16 shows that 
the set of rotor waves and the set of stator waves are the same set. Thus, we have found all waves 
generated by coupled rotor/stator interaction (equation A-13); they are the same set found in uncoupled 
analyses such as that of Tyler and Sofrin (ref. 1)! Furthermore, comparison of equations A- 14 and A- 15 
shows that the roles of the n and k indices switch when the same wave is viewed in different frames. 
In the rotor frame, k is the cascade index and, in the stator frame, it is the time harmonic index. 

Study of the above equations also shows the origin of frequency scattering, equation A- 14 gives the 
form in the rotor frame of the waves caused by rotor loading. A single loading harmonic n produces 
the waves indicated where k takes on all integer values. That is, one loading harmonic produces waves 
at the same frequency but with many mode orders. Equation A- 15 represents these same waves in the 
stator frame and shows that the waves caused by a single harmonic (n) of rotor loading appear at all 
harmonics ( kBjQ . ) of blade passing frequency in the stator frame. 

Kinematic phases for rotor and stator waves can also be expressed in duct-fixed coordinates via 
equation A-3. This result plus the other results derived above are summarized below. 


Waves due to stator load harmonic n and cascade index k 


$nk ( \>2 +nB l Qt 


kB 2 ^ 2 ) t 


in stator frame 
in rotor frame 
in fixed frame 


(A- 17) 


Waves due to rotor load harmonic n and cascade index k 




+nB 2 Qt in rotor frame 

'fifafyi+kBiClt in stator frame 

+ (A:fi| £4| - nB 2 0 . 2 )t in fixed frame 


In code CUP2D the f} nk ’s are computed in subroutine ALFBET. 


(A- 18) 
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APPENDIX B 

AXIAL WAVENUMBERS AND THE CUTON CRITERION 


The general 2-dimensional waves studied in this report are in the form of complex coefficients 
multiplied by the exponential exppYcut+Py+cw)]. Appendix A showed how the frequency co and 
transverse wavenumber p are related to the kinematics of the problem, that is to the numbers of blades 
and the rotational speeds. The axial wavenumber a was derived by Smith from aerodynamic analysis 
as a function of P and 0). Its form determines whether waves propagate or decay in the axial direction. 
In this appendix, we find non-dimensional forms for the axial wavenumbers that can be used for both rotor 
and stator waves. We also evaluate the cuton criterion and deal with a sign choice for the upstream and 
downstream going pressure waves. 

In appendix A we saw that the kinematic phase, \|/ = py + oor took on different forms in the duct-fixed 
and rotating coordinate systems. The axial wavenumbers, however, must be the same in any of these 
coordinate systems. We choose to evaluate them as related to stator waves in the duct-fixed frame. The 
form for the rotor waves will be found simply by permuting indices using the principles derived in 
appendix A. Consider first the wavenumbers for the pressure waves given by equation 4 of the main text 

„ _ I/(to + VP) +a^(co + Vp) 2 -(a 2 -I/ 2 ) p 2 (B-l) 

a 2 -U 2 

From equation A- 17, the frequency in the duct-fixed system is d)=nBjQ.j-kB 2 Q 2 . We define a non- 
dimensional form of this as 

Q nk = (rtBjQ, -kB 2 Q 2 )R/a = nB l M yI ~kB 2 M y2 (B-2) 

where M yl and M y2 are the rotational Mach numbers of the two blade rows. We further define the 
axial and swirl Mach numbers to be 


M x = Via 
M s = V s la 


(B-3) 


Substituting these plus the mode order definition from equation A- 11 into equation B-l gives results in 
the following form for a lnk and a 2nk for the stator pressure waves 


a \,2nk ~ a S l,2 R = 


2,„2 


MJfk* +M S p„*)±J - M s P„/ - ( 1 - M x )p nfc 


(B-4) 


i -m: 


The sign choice is discussed below. For the rotor, it is easily shown that aft = a 12nk permitting the 
forms for a ]nk and a 7nk defined in equation B-4 to be used for both the rotor and stator. 

Equation B-4 contains the desired information regarding cuton. When the square root is real, the 
wavenumbers (upstream and downstream) are real and the waves propagate undiminished. When the 
square root is imaginary, /out in the exponential has a negative real part and the waves decay. Thus, the 
argument of the square root supplies the cuton criterion 

( Q n* +M .tP«fc) 2> ( 1 - M jc) 2 PL for cuton (B ' 5) 
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This is more conveniently expressed in terms of the cutoff ratio t, , which is defined as follows 


Qik + M s 


(B-6) 


which must be greater than zero for a mode to be cut on. This means that, for a mode with $nk lobes, 
the frequency, that is Q nk , has to be "high enough", to satisfy equation B-6, for propagation. If $ nk is 
positive, then the effect of swirl is to lower the RPM for cuton (assuming M s > 0). For fans designed 
for cutoff, that is with B 2 > B } this will be the case because the k = 1 mode is the first to propagate 
at BPF so that $ nk = -( nBj-kB -, ) = -{B r B 2 ) = BtB } >0. Study of equation A-8 and the following text 
reveals that this condition corresponds to the case where the mode is spinning in a direction opposite to 
the swirl, that is the velocity of the mode relative to the fluid is increased by the swirl. 

The sign choice in equation B-4 is made from physical reasoning which is different for decaying and 
propagating waves. For decay, the requirement is simply that the waves diminish to 0 as x — > ±». For 
progagation, we require that the wavelength (reciprocal of a) in the upstream region be smaller than the 
wavelength in the downstream region. The working forms for the axial wavenumbers are summarized 
below using the definitions 

F = (B?) 

which apply in the region between the rotor leading edge and the stator trailing edge. The wavenumbers 
for the vorticity waves are also given in terms of quantities already defined. 


Propagating pressure waves ( E<() ) 


F>0 : 


a \nk ~ 
a 2nk = 


'[F + sTE] 
1 -M x 


-T_[ f-PF] 

1 -M~ 


(B-8) 


“in*— '-^IF-Pe] 

1-AC , DO , 

F< 0: . * (B-9) 
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Decaying pressure waves ( E >() ) 


a 


ink 


- — [F-i\fE ] 


a 


Ink 


i - m : 


i-m; 


.[F + isfE] 


(B-10) 


Vorticity waves 


a 3 nk ~ 


^ nk + ^s 


M 


x 


(B-l 1) 


As a result of nomenclature developed in appendixes A and B, convenient forms can be used for the 
exponential exp^You+Py+cor)]. For the upstream pressure wave caused by stator loading harmonic of 
order n, the exponential is exp[i(a ]nk x+\\i nk )], where k is the cascade index that takes on all integer 
values. Formulas for \\f nk valid in all coordinate systems of interest are given in equation A-17. For the 
upstream pressure waves caused by rotor loading harmonic of order n, the exponential is 
exptifa^^c+Xj/^j], where xj/^ is given in equation A- 18. For downstream waves, change the 1 in the 
a subscripts to 2’s. 


The axial wavenumbers derived above apply directly to the swirl region between the rotor and stator. 
For the region upstream of the rotor, the notation is continued with a 4nk for the upstream- and 
downstream-going pressure waves. Downstream of the stator, a 8nk and a 9nk apply to the downstream 
going pressure and vorticity waves. To apply the formulas above, set M s =0 and use the axial Mach 
number M x appropriate for the region. Also, because £> nk is based on the speed of sound in region 
b, this must be adjusted for the other regions. The results are 

F a =M xa Q nk a b la a (B12) 

for region a and 

F c =M xc Q nk a b la c (B13) 

£ r =(l-<)P' t -( n »*V"r) 2 


for region c. 

In code CUP2D the a’s are computed in subroutine ALFBET. 
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APPENDIX C 

COMMON FORMS FOR TRANSVERSE VELOCITY COMPONENTS 


In reference 7, Smith derives expressions for v/, v/, and v/, the transverse velocity of 
the upstream-going pressure wave, downsteam-going pressure wave, and vorticity wave caused 
by a unit loading element. The axial velocity and pressure are then computed from the v A s. 
In this appendix, we apply the non-dimensionalizations of this report and derive forms for these 
transverse velocity components that can be applied in the coupled cascade code CUP2D to both 
the rotor waves and the stator waves. 


Starting with pressure waves from the stator, we apply Smith’s equations 23 and 28 to find 


v 


/ . 3 



=f 0? A, 2 sin02 


V 


cos0 2 


(C-l) 


which is written directly for decaying waves. The upper/lower signs apply to the up- 
/downstream-going wave. By reproducing Smith’s derivation, it can be shown that the correct 
interpretation for the square root for propagating waves gives 


where 


1,2 


_JL 

2A- 


t (5 * X 2 sin 02 + 


X, 2 3 cos0 2 


A4 *2 




A 2 = X ,2 + ft 2 + 2A, 2 ftsin0 2 


(C-2) 


(C-3) 


The subscripts on the right side have been added to identify waves caused by stator (blade row 
2) loading waves. M r2 is the relative Mach number at the stator. The reduced frequency and 
transverse wavenumber are 



-Q 2 )c 2 
*2 


(C-4) 


and 


P=pc 2 (C- 5 ) 

In terms of the dimensionless quantities c 2 = c 2 /R and P nk = PR, the reduced frequency can 
be rewritten as 
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where 


X 2 - « A 2 C 2 
_ BfM yJ -M y2 ) 


(C-6) 

(C-7) 


and where M yl =Q. J R/a and M y2 =Q.fR/a are the rotor and "stator" rotational Mach numbers. 
(Recall from appendix A that this analysis applies to counter rotation systems; however, for 
convenience, the downstream blade row is called the stator.) Now, the normalized velocity 
components for the upstream-going pressure waves can be written 


/ 




~fink' nA 2 sin0 2 
V 


\ 

^ nA lfink cosQ 2 

ROOT 


(C-8) 


and for the downstream-going pressure waves 



( 


1 

2A 2 


+ Pnfc +nA 2 sin0 2 

V 


\ 

« A 2p«fc COS0 2 

ROOT 

y 


(C-9) 


Normalization by $ nk is needed for cases when f) nk = 0. In these cases, v//f) nk is still defined 
and is needed to determine the axial component and the pressure from equations 7 and 8. In 
equations C-8 and C-9 the complex square root is defined by 


ROOT = 



for propagation 
for decay 


(C-10) 


and 


E 2 = Ak-M 2 r2 A- 


(C-ll) 


A 2 = -z| = n 2 A 2 + title +2n A 2 $nk sin 0 2 


For the vorticity wave. Smith’s equations 24 and 28 give 


(C- 12) 
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(C- 13) 



The above derivation applies to waves from the stator caused by loading harmonics of order n. 
For rotor waves associated with rotor loading harmonic n , the derivation proceeds along the 
same lines but with subscripts 2 (on the right hand sides of the equations above) replaced by l’s 
and (3 nk replaced by P kn . The result is that common forms can be writted for the Vs as 
follows. 


/ 



1 

2A 


-n Asin0 


n A (3 cos 0 
ROOT 


(C-15) 



ROOT- for propagation ( E <0 ) (C-20) 

-i\[E for decay (E> 0) 

Interpretation of the above notation for waves caused by either rotor and stator loading at 
harmonic n is given in the table below. 
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STATOR 


ROTOR 


a: 

a 3nk 

a 3kn 

P: 

$nk. 

$kn 

6: 

0 2 

0 1 

A: 

B { (M yI -M y2 ) 

B 2 (M yl -M y2 ) 


M r2 

Mr, 


In code CUP2D, V ]t V 2 , and V 3 are computed in subroutine GETVS using equations C-15, C- 
16, and C-17. 
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APPENDIX D 

REFLECTION AND TRANSMISSION OF WAVES AT ACTUATOR DISKS 


In this report the effects of mean blade and vane loading are to modify the mean flow 
environment in which the acoustic waves propagate. Thus, the rotor turns its axial inflow in the 
swirl direction and the stator straightens the flow. Blade loading is considered the superposition 
of a mean loading component and the unsteady loading. Unsteady loading is distributed over the 
blade chords and produces small flow perturbations; however, the steady loading can be large and 
affects the wave propagation properties. Since the steady loading by itself does not make noise 
with subsonic rotors, we consider only its effect on wave propagation. For simplicity, the mean 
loading is concentrated in two actuator disks; one at the rotor leading edge and the other at the 
stator trailing edge. Smith’s theory can be applied in any region where the mean flow is 
constant. 

To treat the environment sketched in figure 7, we need to establish the boundary conditions 
(or jump conditions) at the disks. This done in this appendix by deriving reflection and 
transmission coefficients for waves incident on the two disks. As shown in the figure, the 
upstream-going pressure wave, v 1 can transmit a pressure wave v 4 , and can reflect pressure 
and vorticity waves v 2 and v 3 . Similar comments apply at the exit actuator disk. Boundary 
conditions are derived by applying conservation of mass, axial momentum, and transverse 
momentum at the disks for each type of wave on a mode by mode basis. These lead to reflection 
and transmission coefficients as defined in equations 16, 17, and 18. 

Horlock (ref. 12) in his book on actuator disk theory gives the (unlinearized) 2 dimensional 
conservation equations as follows 


Pi u \ - ?2 u 2 

continuity 

(D-l) 

-L cosoc m +p 1 « 1 v 1 =p 2 «2 v 2 

transverse momentum 

(D-2) 

L sina m +p, +p,« 1 2 = p 2 +p 2 u 2 

axial momentum 

(D-3) 


where the definitions for the p’s, u’s, and v’s are standard. The subscripts 1 and 2 are 
Horlock’s and apply to conditions on the two sides of the disk in the above equations only. They 
should not be mistaken for 1 and 2 subscripts elsewhere in this report. The terms with the L’s 
represent axial and transverse mean loading on the disk. These terms will be seen to drop out 
of the analysis via the linearization that separates steady and unsteady effects. In the following, 
the procedure is to linearize the continuity and momentum and equations, express the fluctuating 
quantities in terms of pressure and velocity components, and then apply the notation from Smith’s 
report to represent each wave by its transverse velocity component. As indicated by Figure 7, 
there are three interactions to be considered, each involving an input wave and three resulting 
waves. For instance, at the upstream interface, the input wave v } produces the reflected waves 
v 2 and v 3 and the transmitted wave v 4 . Thus, for these four waves, linearized versions of the 
three equations above can be used to express the resulting three waves in terms of the input 
wave. 
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Upstream-going Pressure Waves 

In figure 7 the flow field is divided into regions a, b, and c. In this section we consider the 
jump condition at the interface between region a and region b starting with the continuity 
equation. 

Continuity Equation 

Equation D- 1 can be expressed in terms of the mean and fluctuating components as follows. 

(D ' 4) 

from which the first order unsteady terms give 

(D ' 5) 

The mean flow quantities p a , p b , U a , and U b are considered known constants for the analysis. 
In order to deal with pressure rather than density for the fluctuating quantities, we apply p = pa 2 
with the result 


_ V a _ _ U b _ 

Pa u a + ~ Pa = Pb u b + ~ Pb 


(D-6) 


a. 


a i 


In the upstream region, only the upstream-going wave exists. For axial velocity and pressure, 
this can be expressed in Smith’s notation as follows. 


— < a 4 x 
u a= u 4 e 

— ia 4 x 

Pa s P4<' 


(D-7) 


where the remainder of the exponential, exp^Ycor+pyj], is implied and is the same throughout the 
field for each mode/frequency. Between the disks, all 3 waves exist as expressed by 


— ia,x — ia?x — io.>x 

«£ = «1 e + u 2 e + u 7, e 

— / a,x — ia 7 x 

Pb=P\ e + Pl e 


(D-8) 


Substitution of equations D-7 and D-8 into equation D-5 gives the relation between the unsteady 
pressures and the unsteady axial velocities 

_ U a - ___ U b - - 

Pfl«4 + —P4 = P^ M 1 +u 2 +u ?) + — (p 1 + Pl) (D-9) 

a a a b 

at the inlet interface (x=0). Now equations 7 and 8 from the main text can be applied to express 
this in terms of complex amplitudes of the transverse velocity components as follows. 
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(D- 10) 


p„[a 4 - i^(m-{/ a a 4 )] £ = ?„[«, - Zl (m^aj * V s p)] 

a a P a b _ 

*P b \o-2 - ^f» + U b a 2 + V s V)}?l 

ft, P 

a 


This can be non-dimensionalized by multiplying numerators and denominators by the effective 


radius R. The result is 

c 4 v 4 = c, vc 2 v 2+ c 3 v 3 

(D-il) 

where we have defined 

''i -*|/(»»Prt) 

Vj = Vj/(a b P n4 ) 

(D-12) 


V, • v,/(a 6 a w ) 






Ci =p b a b [(l -MZ)a lnk -M b (Q nk +M s $ nk )] 

C 2 =Pb a bU [ ~ M b a 2 nk~ M b( Q nk +M s^ 

C 3 =Pb a b[~$nk) 

C 4 = P»««[(l-M 4 )a 44 t -«»—«»*] 

a a 

We define reflection and transmission coefficients as in equation 16 of the main text 

V 2 -R n Vj 

Vj=R 13^1 

These allow equation D- 1 1 to be written in the following 2 forms 

C 4 V r 4 =(C 1 + C 2 /? 12 + C 3 /? 13 )V 1 


(D-14) 


(D-15) 
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and 


(C 4 r 14 -C 1 -C 3 y? 13 )v 1 =c 2 v 2 (D-i6) 

which will be used later when the momentum equations have been applied. 

Axial Momentum Equation 

This section proceeds along the same lines as the derivation above from the continuity 
equation. Equation D-2, in terms of steady and fluctuating components becomes 

L sin «,„+(/■„ *p a ) +(p„ +uf = ( P b +p b ) *( P() *p 4 ) (U b *U b ) 2 (D- 17) 

After expanding this and eliminating the unsteady densities with p = p / a , 
unsteady equation is found to be 

(1 *M 1 a )p u *2p a U a ii a =(1 *M^p b -2p b U b a b 


the first order 
(D-18) 


As anticipated above, L has dropped out and need not be specified. (The effects of L are seen, 
however, through the mean flow parameters.) The same procedure followed above in arriving 
at equation D- 1 1 lead to 

f 4 v 4 = f, vf 2 v 2 *f 3 v, (D-19) 

F 1 p b a b [(l +M 2 b )(Cl nk +M s P„ t ) - < I -M 2 )M b aj 
F,*p b a„l2M b V nk l < D - 20 > 

-PAld M a o. ink ] 

U-h 


The reflection and transmission coefficients of equation 16 produce the final results of this 
section. 


F^-^F^F^ 3 ) V, (D-21) 

and 

^FJu-F.-F^V^F^ (D-22) 

These will be used to find the forms for R n , R ]3 , and t I4 at the end of the next section. 


Transverse Momentum Equation 

In terms of steady and fluctuating parts, equation D-3 becomes 
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-L cosa m +(p a +p fl )(£/ fl + M a )(V a +v a ) = (p ft +p 6 )(i/ 6 +u b )(V b +v b ) 


(D-23) 


Note that in the inlet region the transverse velocity V a = 0 and in the swirl region V b = V s . 
Then, by the same procedures used above, the first order unsteady equation is found to be 

Pa U aV a = Pb U bVb + Pb V s a b +M b M s Pb (D ' 24) 


With the same definitions used above, this becomes 

g 4 v 4 =g,vg 2 vg 3 Vj 


(D-25) 


with 


Gj -p b a b [-M b M s (£lnk +M s$nk) +M b$nk + tt ~ M b M s a \nk^ 

G 2 = Pb a bl~ M b M s( Q nk +M Al} +M bV > nk + ( [ - M b > > M s a 2nlJ 
G 3 ~ Pb a bf- M b a 3nk~ M s$nk} 

G 4 = P fl a a[ — M aW 

a b 


In terms of the reflection and transmission coefficients 

G 4^4 = ( G 1 +G 2 R 12 +G l R \?) V l 
and 

(G 4 r 14 -G 1 -G 3 /? 13 )V 1 =G 2 V 2 


(D-27) 


(D-28) 


Coefficients for Inlet Disk 

The preceeding equations have been derived in a form convenient for the elimination of the 

V"s. For instance, division of equation D-15 by equation D-21 eliminates Vj and V 4 . 
Continuing this technique leads to equations that can be solved for the inlet transmission and 
reflection coefficients with the following results. 


£ 2 £ 3 Eq £5 

^1 ^5 ~ ^2^4 
£ i £ 3 - EqE 4 

£jE 5 -£ 2 £4 
£?£ 8 -£ 6 £ 9 

£ 4 £ 9 -£j £ 7 


(D-29) 


where the £ coefficients are defined in terms of the previously defined C’s, Fs, and G’s as 
follows. 
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Eo-c 4 f,-c,f 4 

E { = C 4 F 2 -C 2 F 4 
£ 2 =C 3 F 4 -C 4 F 3 
£ 3 =C 4 G! -0^4 

F 4 = C 4 G 2 - C 2 G 4 (D-30) 

£ 5 =C 3 G 4 -C 4 G 3 

E 6 =C 2 G 1 _C 1 G 2 
E-j =C 3 G 2 ~C 2 G 3 

^8 = G 2 ^l ~ G \ F 2 

E 9 = C 3 F 2 -C 2 £ 3 


Downstream-going Pressure Waves at Exit 

For waves reflected at the stator exit actuator disk, the derivations are similar to that for the 
inlet. The major difference is that exponentials containing x do not vanish as they did at x=0. 
For later convenience, the exponentials are absorbed into the reflection and transmission 
coefficients, as will be seen below. As indicated by figure 7 for x=x e , a downstream-going 
pressure wave v 2 reflects an upstream-going pressure wave v ; and transmits a pressure wave 
v 8 and a vorticity wave v 9 . Coefficients for these processes are derived next. 

Continuity Equation 

By analogy with equation D-6, the linearized continuity equation is 

u b u c 

Pb“b + —Pb = Pc a c + — Pc (°- 31 ) 

a u a. 


Application of equations D-8 and similar equations for region c at x=x e leads to 

— U u — in.r — U h — in.Y — _ 


(P b 


- U b~^ ia,x e , - u b -x iu 2 x f , - u c-, ia 8 x e ~\, i a <?4D-32) 

u^-lpje 1 e +(p b u 2 +—p 2 )e " = (p r w 8 +__p 8 )e +(p c Ug)e 




a. 


As in the previous section, equations 7 and 8 from the main text are used to express this result 
entirely in terms of the transverse velocities. Then, with the same normalizations used above, 
we find 


Cj e * e +C 2 V 2 e l a2nk Xf 


i a 


= C 8 Vg e ' e + C Q V a? 


l ®-9nk x e 


' 9 Y 9 


(D-33) 


where C ] and C 2 have already been defined, 
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(D-34) 


V 8 =v 8 /(a c $ nk ) 
V 9 = vg/(a c a 9nli ) 


and 


C 8 = p A Kl-M f 2 )a g ,,»-^— 

a c 


(D-35) 


^9 P c a c^ ^ 

The transmission and reflection coefficients are defined with the axial exponential included 


1/ o i a lnk x e _ n \/ 0 1 a 2nk x e 

V _- R 2\ V 2 e 

i7 ia 8nk x e t i7 ia 2nk x e 
V 8 e anK e = f 2 g V 2 £ f 

V 9 e ia ™ x * = T 29 V 2 e ia2 ^ 


(D-36) 


When these definitions are applied to equation D-33, the following two equations can be derived 

(D-37) 


(C,« 2 | *c 2 -c 9 r 29 ) = c 8 v'"**'' 


c, v, e '“'-a = (c 8 r 28 * c, r 29 -c 2 ) v 




for use below. 


Axial Momentum Equation 

Proceeding along lines similar to the above derivations for the axial momentum equation, we 
define 

Fg =p c <J c [(l -»l) — **<•«»] 

(lb (D-38) 

F 9 =Pca c [2 — M^ nk ] 
a b 

and use the coefficients defined in equation D-35 to establish 

F, V, , '“'" A = (F 8 7jg-F 2 -F,r 2 ,) V 2 e (D . 39) 

(F|R 2 | *F 2 -F,T 29 ) V 2 <- io » lr ' = Fg Vgc 

which also will be used shortly. 
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Transverse Momentum Equation 


Here we define 


Gg = p c a c [ — - M c P,^] 

a b (D-40) 

a c 

Gg = p c a c [ — M c 

a b 


From which it follows that 

(G,* 21 ♦ C 2 -C 9 T 29 ) V 2 c ia2 '*' t ' = G 8 V 8 e ,as "*** 
G, V, e io '”‘ x * = (C 8 r 28 *G 9 T 29 -C 2 ) v 2 * ,a! "' x * 


(D-41) 


Now equations D-37, D-39, and D-41 can be manipulated to find the desired coefficients as 
follows 


^21 “ 


E \5 E \6 


E [4 E ll 


T 2S = 


T 29 ~ 


E [2 E \5 E \0 E \1 
E (> E U ~ E % E l?> 

^10^13 “^11^12 
£ 6 £ 10 ~ £ 8 £ 12 
£10^13 -E n E n 


(D-42) 


where 

*10 = ^1 ^8 "^ 8^1 
£jl = CgFj -Cj Fg 
£12 = Cj G 8 -CgGj 

^13 = ^9^1 "^jGg (D-43) 

^14 = ^S E 2 ~^2 E i 
Cl 5 = CgFg “CgFg 
Ci6 = CgG 2 -C 2 G g 

Cj7 = CgGg -CgGg 


Vorticity Waves at Exit 

For this interaction, the derivation is the same as for the pressure wave at the inlet except the 
input wave is different. We quote only the results. The reflection and transmission coefficients 
are 
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^31 = 


^38 ' 




E l 5 E l9~ E n E \S 
E n E \s~ E io E n 

E n E 2l ~ E 13 E 20 
^ 11 ^ 12 -^ 10^13 
£10^19 -£12^18 
E 10 E [1~ E \2 E \5 


(D-44) 


where 


^18 = ^ 8^3 ~ E 2 E H 
E \9 = ^ 8^3 "^ 3^8 
E 20 = E \ E 2 E l 
E 2\ ~ ^ 1^3 "^ 3^1 


(D-45) 


and where the C’s, F' s, and G’s and the rest of the E’s have already been defined. 


This completes derivation of the reflection and transmission coefficients for waves impinging 
on the actuator disks. The coefficients are given in equations D-29, D-42 and D-44 in terms of 
coefficients defined by equations D-13, D-20, D-26, D-35, D-38, D-40, D-30, and D-43. These 
are all computed in subroutine RTCOEF in the code CUP2D. 
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APPENDIX E 

FORMULAS FOR VISCOUS WAKES 


In the main text, the method of inverting the matrix equations for the coupled loading was 
described in conjunction with equation 56. This equation requires a vector of upwash velocities 
WWS(n,i) at the stator caused by viscous wakes from the rotor, n denotes the harmonic order 
and i denotes the control point on the stator. This appendix derives the working form of these 
formulas for use in code CUP2D based on the viscous wake correlation of Silverstein, Katzoff, 
and Bullivant (ref. 13). 

Figure E-l shows the variables for the derivation. The objective is to start with a spatial 
description of the wake velocity defect fixed in the rotor frame, convert this to a waveform of 
upwash in the stator frame via coordinate transformation, and then Fourier analyze this to find 
the desired harmonics. From the sketch with the airfoil and wake in Figure E- 1 , the form of the 
defect is given by 

— =exp[-it (y/Y ) 2 ] 
w c 

where the maximum defect is given by 

_ 1.21 

w i Zj -0.70 

and the wake width parameter by 


(E-l) 


(E-2) 


— = 0.68 ^(zj -0.85) (E-3) 

c i 

The last two formulas give the decay in magnitude and spread in width as a function of 
chordwise distance measured from the leading edge, t is normalized chordwise distance z l /c 1 
and C D is the standard isolated airfoil drag coefficient defined by C D = (drag/unit 
span)/(’/ 2 pI/ 2 *chord). 


Wake decay and spread are to be specified at the point where the wake is intercepted on the 
stator. This can be found by relating the chordwise positions to the x direction as follows 

Zj cos 0j = x s + Z 2 cos 0 2 (E-4) 

In terms of normalized quantities, this leads to 


x s + z 2 c 2 cos 0 2 
c j cos 0 , 


(E-5) 


To obtain the waveform in the stator frame, we note that the relative motion of the two bla 
rows can be specified in terms of y 7 and y 2 , and then we locate the position in the wake 
in terms of y l for a fixed x as follows 
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v \a> 


And since 


y ( =xtan0j + y/cos0j 


(E-6) 


I y 2 =yi ~y\ (E-7) 

it follows that 


y =y 2 cos0j -cos0j QRt -.xsinfij (E-8) 

By expressing x and y 2 in terms of z 2 and substituting into equation E-l, the wake waveform 
in the stator frame can now be written 


w ’ v (t) = w c fit) (E-9) 

where 


(E-10) 


If we recognize that this one pulse in a periodic waveform of period In/Bfo, it is a 
straightforward matter to express fit) in a Fourier series of the following form 

(E-io 

n=l 

where, as usual, the real part is implied. The usual Fourier integral formulas lead to 

F n = 2 ^~ zxy[-{KnlQ ) 2 ] g inB i [ ~ 2 2~2 s > n Ce 2 -e 7 ) + J,sin0 ; ]/cos0 / (E-12) 

where 


fit) = exp 
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sin(02-0j)z2 -JCjSinOj -cos0j QRt 


Q=\[k 


5C, 


cos0. 


(E-l 3) 


Now, we have only to find the component of w v normal to the stator chord. This is given by 
w v sin(0 2 -0j). Hence, the desired upwash harmonics can now be written 


w 


VS 


Wo 


Jn 


vt> 


W. 


COS 0O 

sin (8 2 -6,) 


(E- 14) 


where, for convenience, we have used WjAV 2 = cos0 2 /cos0,. When w/Wj and Y/c l are 
evaluated from equations E-2, E-3, and E-12 at values of t 2 corresponding to the control points 
Zi used by Smith, the result is WWS(n,i) as needed for inversion of the matrix equation. 

The above formulas are equivalent to those used by Kaji and Okazaki except that we retain 
the absolute phase. In code CUP2D, the above formulas are coded in subroutine GTWAKE. 


i 
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APPENDIX F 
LIST OF SYMBOLS 


a = local speed of sound 

c = blade chord 

Cj,c 2 = rotor, stator chord 

i - index to count control points on blades and vanes 

j = index to count loading panels on blades and vanes 

k = cascade index 

m = nBj-kB 2 , circumferential mode order 

n = loading harmonic index 

p = disturbance (acoustic) pressure 

u = axial component of disturbance velocity 

v = transverse component of disturbance velocity 

v = Smith’s basic solution for the transverse velocity, see appendix C 

s = blade gap (spacing between blades or vanes in y direction 

t = time 

w r = maximum velocity defect on viscous wake centerline 
w v = velocity profile in viscous wake definition 

x - axial coordinate 

x e = axial distance from leading edge of rotor to trailing edge of stator 

x s = axial distance from leading edge of rotor to leading edge of stator 

y = transverse coordinate 

z o> z oj ~ source coordinate measured in chordwise direction (at angle 0 to axial direction) 


Bi = 

B 2 = 

BPF = 
L 

LR, LS = 

M 

M r 

M s = 

M yl = 
My2 

»h = 

%l : 

R 

R 12 ,e tc = 
SC 

T I4 ,etc = 

U 

V 

Vs = 

w 

WR,WS= 


number of blades on upsteam blade row (which is usually called the rotor) 

number of blades on downsteam blade row (which is usually called the stator) 

blade passing frequency 

discretized loading source 

loading elements on rotor, stator 

Mach number 

Mach number relative to blade, W/a 

Mach number of the mean flow in the swirl (y ) direction 

rotational Mach number of upstream blade row (normally the rotor) 

rotational Mach number of downstream blade row (normally =0 for the stator) 

number of loading harmonics included in coupling 

number of panels on rotor (same as number of panels on stator) 

sound power level in dB, see equation 1 14 

effective radius of annulus, used to convert transverse linear coordinate to angle, <) )=y/R 
modal reflection coefficient, see equation 16 and appendix D 
spacing to chord ratio, s/c 

modal transmission coefficient, see equation 16 and appendix D 
axial component of mean (background) flow 

transverse component of mean (background) flow, relative to the blades 
transverse component of the mean swirl flow, absolute frame 
[U 2 +V 2 ]' /2 = mean flow speed 
upwash velocity (normalized) on rotor, stator 


Note, the analysis of this report applies to counter rotation stages; however, for convenience, 
the upstream blade row is called the rotor and the downstream blade row is called the stator. 
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a = 

p 

$nk 

0 

(0 = 

Pa’Pb’Pc= 
Po 

a = 

¥ 

¥„* = 


axial wavenumber 
transverse wavenumber 
- ( nBj-kBj ) 

cascade angle (determined by U and V) 

angular frequency 

mean density in regions A, B, C 

local mean density 

interblade phase angle 

angle in cylindrical coordinates 

kinematic phase, see equation 1 1 

see equations A- 17 and A- 18 (has different forms in different coordinate systems) 


r = vorticity (loading source) 

Q = angular speed of rotor 

Q nk = non-dimensional frequency, see equation B-2 


Subscripts 

1 denotes upstream-going pressure wave in region B, also denotes rotor 

2 denotes downstream-going pressure wave in region B, also denotes stator 

3 denotes vorticity pressure wave in region B 

4 denotes upstream-going pressure wave in region A 

5 denotes downstream-going pressure wave in region A 

6 denotes vorticity pressure wave in region A 

7 denotes upstream-going pressure wave in region C 

8 denotes downstream-going pressure wave in region C 

9 denotes vorticity pressure wave in region C 

a denotes condition for region A (upstream of rotor leading edge) 
b denotes condition for region B (between rotor leading edge and stator trailing edge) 

c denotes condition for region C (downstream of stator trailing edge); also used in w r , the 

maximum velocity defect on the wake centerline. 
i = index to count control points on blades and vanes; also, imaginary unit 
j = index to count loading panels on blades and vanes 

k = cascade index 

n = loading harmonic index 

r denotes rotor 

rr denotes effect of rotor on rotor 

rs denotes effect of rotor on stator 

sr denotes effect of stator on rotor 

ss denotes effect of stator on stator 

s denotes tangential (swirl) component, also denotes stator 
x denotes axial component 

R, I denote real and imaginary parts of a comlex quantity 

overbar denotes normalization by effective radius, R for x’s and c’s 
( ) denotes complex conjugate 
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Figure 1 . Three-dimensional rotor and stator geometry (ref. 3). 
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ROTOR 


STATOR 



Figure 2. Simple rotor/stator interaction model. Dashed blades indicate that their role is only 
to generate wakes. They are transparent to acoustic waves. 



FAN TIP SPEED, FT/SEC 


Figure 3. Aft radiated sound power levels from engine designed for BPF cut-off. (Topol, et al . , 
ref. 5). 
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Figure 4. Two-dimensional geometry of Smith’s theory (ref. 7). 
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Figure 5. Geometry for coupled cascade analysis shown at t= 0. Either blade row or both blade 
rows can rotate. Note that 9, is normally negative. Also, if the downstream blade 
row rotates, Q 2 would be negative. 
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ROTOR STATOR 

INLET EXIT 

DISK DISK 


Figure 6. Direct waves ( KS1 , KS2, KS3 ) from load element LS( nj) and waves reflected ( VS 1 , 
VS2, VS3 ) from actuator disks. Note that the entire wave system is driven by LS(nj). 
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Rotor 
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disk 


Stator 

exit 

disk 





x 


Figure 7. Three types of reflection/transmission included in the method. Vj, v 2 , v 4> are the 
transverse velocity in the pressure waves (denoted by solid arrows). The dashed 
arrows denote vorticity waves with transverse velocity v ? and v 9 . 
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Figure 8. Comparison with Kaji and Okazaki theory for IGV/rotor interaction. Gap/chord ratio 
= 1.0 for both blade rows. Also, Bj=B 2 Drag coefficient = 0.02. 
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Figure 9. Computed results of Kousen and Verdon (ref. 1 1 ) using LINFLO for stator excited 
by vorticity wave. Plots shows the resulting pressure wave, which is cut on upstream 
and cut off downstream. 
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Aft PWL Aft PWL 



0.65 0.70 0.7S 0.60 0.85 0.90 0.95 1.00 

Rotor Rotational Mach Number - M yJ 

a) Effect of coupling for one harmonic only 
without swirl. 



Rotor Rotational Mach Number - M y i 


b) Effect of coupling (and frequency scattering) 
without swirl. 



0.65 0.70 0.75 0.80 0.85 0.90 0.95 1 .00 

Rotor Rotational Mach Number - M y j 

c) Effect of swirl for one harmonic with coupling. 
(Note: Major effect is that cascade angles 
change with swirl angle.) 



0.85 0.70 0.75 0.80 0.85 0.90 0.95 1.00 

Rotor Rotational Mach Number - M y i 
d) Effect of coupling with swirl. (All effects included.) 


Figure 10. Computed results showing the effects of unsteady coupling, frequency scattering, and 
swirl on the harmonic power levels propagated downstream. B } - 38, B 2 =72 y A^=30, 
BVGAP= 2.0. Excitation by wake fundamental harmonic only. PWL defined in 
equation 114. 
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REFLECTION COEFFICIENTS AT INLET INTERFACE 



ROTATIONAL MACH NUMBER - M ^ 


Figure 11. Absolute value of reflection coefficients for n,k mode. a,b,c denote cuton speeds 
for the mode in the upstream, interstage, and downstream regions. 
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OF BPF 

• TRAPPED BPF 
WAVE AMPLIFIES 

• STRONG FEED 
TO 2 & 3 X BPF 



2 & 3 X BPF 



Figure 12. Physical model of mode trapping mechanism. 
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9 0 



Rotor Rotational Mach Number 


Figure 14. Comparison of harmonic power levels upstream and downstream. fi 7 = 38, B 2 =72, 
BVGAP=2.0, N p =30, N h =3. Arrows along abcissa indicate cuton speeds for relevant 
modes. Notation is n, k r where n, k indicates mode order and r is b or c 
denoting the region (inter-row, or downstream) as shown in figure 7. 
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0.65 0.70 0.75 0.80 0.85 0.90 0.95 1 

Rotor Rotational Mach Number 

Figure 15. Predictions with and without wake upper harmonics set to 0. Bj= 38, B 2 =12, N = 30. 
N h = 3, BVGAP= 2.0 













Rotor Rotational Mach Number 

Figure 16. Effect of changing vane numbers at constant solidity. Bj= 38, N p =30, N h = 3, 
BVGAP= 2.0. Notation on abcissas for mode cuton points same as in figure 14. The 
*’s indicate trapped modes. 


75 










BVGAP= 1.0 


BPF 



Figure 17. Effect of changing rotor/stator spacing. S i =38, B 2 =12,N p =30, N h =3. BVGAP=axial 
distance from rotor trailing edge to stator leading edge in rotor chords. 
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